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ABSTRACT 


This thesis presents the formulation of structural dynamic equations of motion 
of a general bearingless rotor blade configuration with swept tip The equations 
of motion have been derived using Hamilton’s principle and Rayleigh-Ritz finite 
element discretisation of the beam First the formulation is validated by comparing 
the results of the present analysis for a straight bearingless rotor blade with those 
available in literature 

Then a systematic structural dynamic analysis has been performed to identify 
the effect of control system, location of elastomer and flexbeam-blade-torque tube 
junction from root The effect of tip sweep angles on the natural frequencies and 
mode shapes of the beanngless rotor blade, has also been analysed 

Quantitative measures of lag-torsion and fiap-torsion couplings have been defined 
based on the contribution of torsional deformation in the in-plane bending and 
out-of-plane bending modes respectively The effect of tip sweep angles on these 
couplings have been identified. 
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Chapter 1 


INTRODUCTION 


The dynamic characteristics of rotor blades play significant role in the design of 
helicopter The main rotor system has to provide lift, thrust and control forces 
Therefore, any analytical study pertaining to the dynamics of helicopters requires 
the development of suitable mathematical models for 

Rotor blade 

Fuselage 

and Rotor-fuselage interface 

Rotor blade model essentially consists of the development of structural, iner- 
tial and aerodynamic operators associated with its motion The fuselage model is 
represented by an idealised structural model of a 3D structure The rotor-fuselage 
interface model must represent both the geometry of the interface as well as the aero- 
dynamic interaction in an appropriate manner The fundamental tool, necessary for 
the analysis and design of helicopters, is the rotor blade dynamic model 

During operation, the blades experience large bending and centrifugal loads In 
order to relieve the rotor bending moments experienced by the blades, early rotor 
blades were provided with flap (out-of-plane bending) and lag (m-plane bending) 
hinges at the root of the blade In addition, a pitch control bearing was provided 
to control rotors Such rotor systems are usually referred to as articulated rotors 
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A schematic diagram of an articulated rotor system is shown m Fig 1 1 The 
large number of moving parts leads to a mechanically complex rotor hub system 
accompanied by the associated wear out problem requiring frequent maintenance 
and replacement of parts With advancement in technology, increasing emphasis 
has been placed on the development of hingeless rotor systems The construction 
of these rotors are relatively simple, because of the absence of flap and lag hinges, 
but a pitch bearing is still provided for blade pitch control A schematic diagram of 
hmgeless rotor system is shown in Fig 12 

The development of a bearingless rotor system aims to eliminate both the pitch 
bearing as well as the external damper by incorporating a highly specialised elas- 
tomer In this rotor system, the blade is attached to the hub through a flexible 
structural element called flexbeam The flexbeam is designed to provide the re- 
quired stiffness in flap and lag deformations of the blade, but it is highly flexible in 
torsion Surrounding the flexbeam, there is a stiff cuff denoted as torque tube which 
IS attached to the blade-flexbeam junction at the outboard end and to a pitch link at 
the inboard end as shown in Fig 1 3 The pitch control of the blade is achieved by 
rotating the torque tube through up/down movement of the point P which in turn 
twists the flexbeam An elastomer is placed between the torque tube and flexbeam 
to provide adequate lag damping It also serves the purpose of a spacer between 
the torque tube and the flexbeam Though a bearmgless rotor is mechanically sim- 
ple, its dynamic analysis becomes very complicated due to presence of elastomer, 
multiple load path and existence of a kinematic constraint at the pitch link 

Eversmce the first successful flight of a helicopter in 1942, there has been a 
continued research effort to develop mechanically simple mam rotors with acceptable 
performance and stability characteristics The development of bearmgless rotors 
aims to achieve this goal Recently bearmgless mam rotors have been developed by 
Boeing-Sikorsky company in US and by Eurocopter in Germany These vehicles are 
currently undergoing flight tests The details of the development of rotor blades can 
be found in Ref [1] 
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In order to develop an indigenous rotor blade dynamic model, m Ref [2], Punit 
Kumar Gupta formulated a finite element model of a hmgeless rotor blade under- 
going flap, lag, torsion and axial deformations The rotor blade was treated as a 
straight uniform beam Subsequently, in Ref [3], Venu Gopal extended the model 
by including non-uniform properties of the blade and swept tip These rotor models 
included all the complex geometric parameters such as torque offset, blade root off- 
set, pre-cone angle, pre-droop angle, pre-sweep angle, pretwist, tip sweep angle and 
tip anhedral angle 

Recently, an idealised and highly simplified beanngless rotor model was develope- 
d by Pohit [1], wherein the elastomer was represented as a combination of nonlinear 
spring and damper In a preliminary analysis, the elastomer was treated as a linear 
spring Then, using Power series expansion and Rayleigh-Ritz method the uncou- 
pled natural frequencies and mode shapes of a rotating blade in isolated lag and flap 
modes were evaluated 


1.1 STRUCTURAL MODELLING 

Since the rotor blades are long, slender beams, they will undergo moderate defor- 
mation A non-linear strain-displacement model is used to describe the coupling 
effects between axial, bending and torsion deformations Generally, the strains are 
assumed to be small in comparison to unity Such an assumption is consistent with 
the design requirement based on fatigue life consideration which states the rotor 
blades must be designed to have an operating strain level well below the elastic 
limit of the blade material 

During the derivation of equations of motion of the rotor blade with swept-tip, 
the inertia and structural operators generate a large number of higher order non- 
linear terms So one has to use a consistent ordering scheme to avoid complexity in 
formulation by assigning orders of magnitude to all the non-dimensional parameters 
of the problem, interms of blade bending slope (which is assumed to be of the order 

3 


I 



of e ) A third order approximation is used in the formulation i e , the terms of order 
are neglected in comparison to terms of order 1, i e , 

1 + 0{e^) « 1 


1.2 OBJECTIVES 

Objectives of present study are 

(1) to develop a most general beam type finite element model for bearingless rotor 
blade, including all complex geometric parameters, such as torque offset, blade root 
offset, pre-cone angle, pre-droop angle, pre-sweep angle, pretwist, tip sweep angle, 
tip anhedral angle, elastomer, torque tube and control spring, 

(2) to validate the model by comparing the results of the study with those avail- 
able in literature, 

(3) conduct detailed studies on the dynamic behavior of bearingless rotor blade 
with straight and swept tips to determine the effects of bending-torsion coupling 
due to tip sweep, anhedral angle and control spring 

A schematic diagram of a bearingless rotor blade with tip sweep is shown in Fig 
1 4 The present study on the dynamic analysis of bearingless rotor blade, is an 
extension of the work performed by Venn Gopal on hingeless rotor blade [3] In this 
study, the nonlinear elastomer is idealised as a combination of three linear springs, 
one each for flap, lag and torsional modes respectively 
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(b) Actual rotor hub (from Ref. [l]). 


Fi(;uic 1 3' ScluMnalic thagisun of Bcminglcss 


Rotor without swept tips 




Section A-A 



Figure 1 4 Schematic diagram of Bearingless Rotor with swept tips 
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Chapter 2 


ROTOR BLADE MODEL AND 
ASSUMPTIONS 

2.1 ROTOR BLADE MODEL 

Dynamics of a bearingless rotor blade is a complex phenomenon In the case of 
articulated and hingeless rotor systems, the air/ground resonance is usually avoided 
by providing external mechanical dampers m the lead-lag mode of rotor system 
Whereas in case of bearmgless rotors, a specialised elastomer is used as damper 

The geometric parameters describing the configuration of the bearmgless rotor 
system is shown in Fig 2 1 The parameter a represents the torque offset which 
IS the distance from the center of rotation (hub center) to the reference axis of the 
blade ei and 63 refer to the blade root offset distance from hub represents 
precone angle defining the orientation of the blade pitch control axis with respect to 
the hub plane Ps and pd represent the presweep and predroop angles respectively, 
representing the inclination of the blade reference axis with respect to the pitch 
control axis The bearmgless rotor can be divided into three parts, namely, (1) 
flexbeam, (2) blade and ( 3 ) torque tube The flexbeam-blade-torque tube junction 
IS at a distance from root and the free end of the torque tube is at a distance ^2 
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from the root An elastomer is placed in between flexbeam and torque tube at a 
distance of Lgi from the root The blade consists of a straight portion and a swept 
tip whose orientation relative to the straight portion is described by sweep angle As 
(aft-sweep is positive) and anhedral angle Aa (positive upwards) A control spring 
(Kc) IS attached to the torque tube at a distance Lc from root through a pitch link 
with an offset Uc 

2.2 BASIC ASSUMPTIONS 

In the formulation of the dynamic model of the rotor blade with swept tip, several 
assumptions have been made which are given below 

1 The blade is treated as an elastic beam 

2 The blade is modeled by beam type finite elements along the elastic axis of 
the beam 

3 The speed of rotation (0) of rotor is constant 

4 The rotor shaft is rigid 

5 The blade undergoes moderate deformation m flap, lag, torsion and axial 
modes 

6 The blade cross-section can have general shape with distinct shear center, 
aerodynamic center and center of mass 

7 The blade is assumed to be made of isotropic material but can have nonuni- 
form properties along the span 

8 The elastomer is assumed as a combination of three linear springs, one each 
for flap, lag and torsion modes respectively 

9 The effect of control system is represented by a linear sprmg(Kc) attached at 
pitch link as shown in Fig 2 1 
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2.3 ORDERING SCHEME 


In the formulation of the equations of motion of a rotor blade with swept tip under- 
going moderate deformation, a large number of higher order terms are generated 
In order to identify and eliminate higher order terms in a consistent manner, an 
ordering scheme is employed This ordering scheme is based on the assumption that 
the bending slopes of the deformed elastic blade are moderate and of order e In the 
derivation of governing equations, it is assumed that terms of order are neglected 
with respect to terms of order 1, i,e , 

1 -h 0(e3) « 0(1) 

The order of magnitude of various non-dimensional parameters governing this 
problem are given below 

ORDER 1 

cos4>k,sin<j)k,Aa,K = 0{1) 

f=0(l) 

hiO = iO = oii) 

ei,6,ac = 0(l) 

ORDER ei 
Ogk = 0(e2) 

ORDER € 

Ql ^ ^ R Ilk. VlX ^ 0(f\ 

Vk,'Wk^(l>,l3p,PdJs = 0(c) 

ORDER cl 
/m,,,,, = 0(c2) 
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^yi — 0{e'^) 


ORDER 

ujt, rnrim, rnCm = 0{e^) 

It IS important to note that ordering scheme is based on physical understanding 
of the behaviour of actual blade configurations Hence care must be exercised in 
deleting higher order terms, based on this ordering scheme 


12 



Blade elastic axis 


Hub axis 
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Hub center 
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(a) FRONT VffiW 



Figure 2 1 Rotor blade with tip sweep and anhedral 
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Chapter 3 


COORDINATE SYSTEMS 


The description of the complex deformation of a rotor blade requires several coor- 
dinate systems The transformation relation between quantities referred m various 
inertial, non-inertial coordinate systems is to be established before deriving the e- 
quations of motion of the rotor blade These coordinate systems are directly taken 
from Ref [3] The relation between two orthogonal systems and Xj,Yj,Zj 

with ext,eyi,ezi and exj,eyj,ezj as unit vectors along the respective axes can be 
written as 


f ^ 



^xi 


Cxj 

By, 

II 

Cyj 

^Zl 

\ J 


\ 


(3 1) 


where the transformation matrix [Ty] can be obtained using the Euler angles re- 
quired to rotate the j-system so as to make it parallel to z-system The coordinate 
systems used in deriving the equation for the rotor model are described below 
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3.1 HUB FIXED INERTIAL SYSTEM-R 


The coordinate system R shown m Fig 3 1, has its origin at center Oh oi the 
unperturbed hub The xr axis is pointing towards the helicopter tail and Zr is 
pointing upwards The unit vectors along the three axes are Cxr, eyR, e^R 


3.2 HUB FIXED MOVING SYSTEM-H 

The coordinate system H shown in Fig 3 2, is a body fixed system with its origin 
fixed at the centre of rotor hub Or Prior to perturbational motion of the hub, 
H-system coincides with R-system If dz,6y and 9x represent the sequential yaw- 
pitch-roll rotations, then the transformation matrix [T^^] can be written as 


[Thr] = 


1 

0 

0 


cos 9y 

0 

— sin 9y 


COS 02 

sm 02 0 

0 

COS^x 

sin^x 


0 

1 

0 


— sm 9g 

sm 02 0 

0 

— sm 9x 

COS^x 


sin^j, 

0 

COS$y 


0 

0 1 


(3 2) 


Since 6x, 9y and dg are assumed to be of order sine and cosine functions can be 
approximated as 


sm 6^9, cos ^ w 1 

Substituting the approximation, [Trr] can be simplified as 


[Trr] = 


1 9g 
9y9g - 9g 1 
9x9z + 9y 9y9z - 9x 


-9y 

9x 

1 


(3 3) 
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3.3 HUB FIXED ROTATING SYSTEM-1 


The coordinate system shown in Fig 3 3, rotates about zih axis with constant 
angular speed Q, of the rotor Its origin is fixed at hub centre Oh This system can 
be obtained by rotating H system by an azimuthal angle ^pk of the blade about 
ziH axis The transformation matrix is given as- 


[Tih] = 


cos -ipk sin ipk 0 
- sin ipk sin ipk 0 
0 0 1 


where, 'tpk is azimuth angle of blade 


'^k — + {k — l)^ and = 


(3 4) 


3.4 ROTATING SYSTEM-2K 

The coordinate system 2K shown m Fig 3 4 is a blade fixed coordinate system 
which rotates with k^^ blade The origin of the 2K system is at the location of the 
blade root A (Fig 3 4), which is at a distance aiyi + eie^ri from the hub centre 
1-system and 2K-system are parallel 


3.5 PRECONED, ROTATING SYSTEM-3K 

The system 3K shown in Fig 3 5a and 3 5b also rotates with blade This system is 
obtained by rotating 2K-system by an angle -j5p(precone angle) about y 2 k axis The 
transformation matrix between 2K and 3K is given as- 
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1 0 

[Tz2]= 0 10 (3 5) 

0 1 _ 

3.6 PREDROOPED, PRESWEPT, PITCHED, 
BLADE-FIXED ROTATING SYSTEM-4K 

The 4K system shown in Fig 3 5a and 3 5b, is blade fixed system with its origin 

at pitch bearing B(or flexbeam-blade-torque tube junction) of the blade It may be 

noted that the pitch axis of the blade is along Cxzk direction and the blade reference 

elastic axis is along the ex 4 k direction While changing the control pitch input of 

the blade, the elastic axis will move on the surface of a cone whose vertex angle is 

described by the angles Ps and Pd as shown in Fig 3 5c 

The 4K-system is obtained by the following steps- 

Step-1 Translating the origin of 3K system by a distance ’62’ along ixzk 

Step-2 Then rotating the system by an angle (presweep angle) about zsk axis 

Step-3 Then rotating the system by an angle (predroop angle) about yzk 

Step-4 Then rotating the system by an angle 61 about xzk 

The transformation matrix is given as- 

1 -{Ps cos 9i + Pd sin 61 ) [Pd cos 0/ - Ps sin Oj) 

[^43]= {PsCosOi + PdSmOj) cos 9 1 sin9j (3 6) 

( - Pd cos 9 1 + Ps sin 9 1 ) - sin 9i cos 9i 
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3.7 UNDEFORMED ELEMENT COORDINATE 


SYSTEM-e 

The e system, shown m Fig 3 6, has its origin at the inboard node of i-th finite 
element The vector ixe, is aligned along the elastic axis of the beam element, while 
the vectors Cyg and e^e are cross-sectional axes For the straight portion of the blade, 
the {ixe, Cyei ^ze) System has the same orientation as {ex 4 k, Cy^k, &z 4 k) system For the 
swept-tip element, the e system is oriented by rotating the Ak system about e^^k by 
sweep angle A^, and then about ty^k by the anhedral angle Aq 
T he transformation matrix between Ak and e system is given as- 


^xe 


* \ 

^xAk 


‘ = [Te4] < 

^y4k * 

^ze 

> 


^zAk 

< / 


For the element in the straight portion of the blade 


(3 7) 



1 0 0 
0 1 0 
0 0 1 


For the swept-tip element 


(3 8) 



COS Ai cos Aq — sin A^ cos A^ sin Aa 
sm Ka cos Aa COS As sin As sin Ao 
sin Aa 0 cos Aa 


(3 9) 


where As is the blade tip sweep angle (positive for backward sweep), and Aq is the 
blade tip anhedral angle (positive upwards) 
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3.8 ROTATING, BLADE-FIXED SYSTEM-5K 


The 5k system, shown in Fig 3 7 is a cross-sectional coordinate system of the 
blade In the undeformed state of the blade, both e and 4k system are parallel 
But the origin of the 5k system is at a distance Xk from the origin of the e system 
During elastic deformation of the blade, the 5k system translates and rotates with 
the cross-section After deformation, the origin of the 5k system, from the origin of 
4K system, is at the location given by 


n— 1 

( ^ ^ (^A: “I” '^k)^xe '^k^ye “t“ "^k^ze (3 10) 

i—l 

The transformation matrix between e and 5k system is obtained following a flap-lag 
sequence of rotation The Euler angles are respectively —pk and Cit corresponding to 
the local slope of the deformed blade in flap and lag directions The transformation 
matrix is given by 



COS Cjt sin Cfc 0 


cos Pk 0 sin Pk 

[r5.1 = 

-sinCfe cosCife 0 


0 1 0 


0 0 1 


- sm Pk 0 cos Pk 


(3 11) 


Since the angle (^f. and Pk are of the order 0(e), the transformation matrix can be 
simplifled by noting that- 


sin 6^0, cos 6 


The Euler angles can be expressed in the terms of the local slope of elastic defor- 
mation of the blade as 


Ck = Vk^^ 
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n— 1 

Note x='Y^ {le)t + Xk 

l~\ 

Substituting the above relations in the transformation matrix [Tse] yields 

1 

■v'k 1 -v,w, (312) 

w'k 0 1 

3.9 COORDINATE SYSTEM-6K 

The 6K system shown in Fig 3 8 represents the cross-sectional coordinate system 
in the deformed configuration of the blade e^, — represents the directions of the 
cross-sectional principal axes 6K-system is obtained by rotating 5K system about 
ex 5 k through the angle {(j)k + 0g), where 6g represents the geometric twist angle of 
the cross-section and (j)k represents the elastic twist 
The transformation matrix relation is given as- 




0 

COS(0G + (t>k) 
sin(^G + h) 


0 

- sin(^G + (t>k) 
cos(6>g -f- <f>k) 



(3 13) 
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Pitch control axis 


Xl, X2k 


(a) Front view 




Figure 3 5 Blade coordinate system 3k and 4k 
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^5k’^z5k 

A 


n A 



Figure 3 7 Rotating blade fixed system - 5k 
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Chapter 4 


KINEMATICS 

During operation, the rotor blade undergoes in-plane bending (lag), out-of-plane 
bending (flap), torsion and axial modes In addition, the hub centre has both 
translational {Rx, Ry, Rz) and rotational {6xi Oy, 9z) perturbational motion The for- 
mulation of inertia operator and aerodynamic operator requires a proper description 
of kinematics of the blade motion In this section, an expression for the absolute 
velocity vector at any arbitrary point ’p’ on the blade is derived 

4.1 POSITION VECTOR OF A POINT 

The position vector of any arbitrary point ’p’ in the finite element of the deformed 

blade with respect to the hub centre Oh, is given by 

n-l 

Tp = aCyi + eiix^k + 62613* + Xl(^e)r6 

t=l 

{Xk “I" Uk^Sxe *1" ^k^Ve '^k^Ze d" ^6jj -b (4 1) 

Transforming all the unit vectors to the 4K-system and neglecting the higher 
order terms, the position vector can be written as 
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n-1 

= ?x4A:^[~^(/^5 cos6j + /3ifSin0/) + ei + 62 + y^Xh)i 

Z=1 

+ cos As cos Aa{xk -\-Uk) + Vkt sin As cos Aq — Wk sin Aq 
+r]{-'^k cos As cos Aa C0 s{6g + <t>k) — w'k cos As cos Aq sin (0(5 + <l)k) 

- sin As cos(0G + <Pk) + cos A^ sin Aa sin(0G + (^fc)} 

+C{UfcCOsA5COsAaSin(0G + ^fc) - WfcCOSAj, COsAaCOs(0G + 0fe) 

+ sin AsSin(0G + <?!'*) + cosA^ sinAaCOs(0G + <l>k)}] 

+ey^J[acos6r + ei(/3sCos07 + ^dSmOj - ^pSmOi) + e2(/0sCOS0/ 

Sin 0/) sin As{xk + Uk) + Vk cos As + r]{-'^k sm A* cos Aa cos(0g + <Pk) 

-w'k sin As cos Aa sin(0G + (l>k) + cos A* cos(0g + <f>k) 

+ sin Ai sin Aa sin(0G + ^fc)} +C{^^fcSinA4COsAoSin(0G + <^jfc) - w^sinAs 
cos Aa sin(0G + ( l > k ) - cos As sin(0G + < f > k )+ sm Ag sin Aa cos(0g + < j > k ) }] 
+ ez ^ J [- a sm 0/ + ei{-fid cos 0/ + 0 s sm - 0 p cos 0/) + e2 {- 0d cos 6 i 
+ 0 S sin 0/) + cos As sm Aa{xk + Uk ) + Vk sm A^ sm Ao + Wk cos Aa 
+' n{Vk sm Aa COs(0G + ( pk ) + fjj'k sm Aa Sin(0G + < i > k ) + cos Aa Sm(0G + ^ k )} 
sm Aa COs(0G + 0 k ) - v'k sm Aa Sin(0G + 0 k ) + cos Aa COS(0G + ^fc)}] 

(4 2) 

In Eq 4 2 all the length quantities are non-dimensionalised with respect to 
the length of the blade I For the sake of convenience, the non-dimensional length 
quantities are referred without an overbar(") The expression 'given m Eq 4 2 is 
general, m the sense that for position vector of any point m the straight portion of 
blade Aa and As are set equal to zero Eq 4 2 can also be written m symbolic form 
as 


Tp — ^[rxex4* + TyGy^^ + Tztz^^ 
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4.2 ANGULAR VELOCITY VECTOR 


The angular velocity vector Co of blade consists of three components They are 

• The constant rotational speed (Cl) of the rotor, 

• the rigid body angular velocity(cnrt 5 id) of the hub due to perturbational rotation 
in roll-pitch-yaw (dx,^y,^z) 

• the angular velocity contribution due to the rate of change of control pitch 
input QOi to the blade 


Ul — "h ^rigid d' 


where 


^rigid — d” ^^y®2/ijt d" ^^z^zik 


(4 3) 


(4 4) 


where ( ) indicates differentiation with respect to the non-dimensional time tp, 

xp = Q,t 


Note that 

Transforming all the unit vectors of Eq 4 3 to the 4K-system and neglecting 
the higher order terms, the angular velocity of blade can be written m symbolic 
form as 

to = d" d” (”^ 

Where 

ujx = [Ox cos xpk d- 9y sin xpk d- 6j {Pp d- Pd cos Oj - Ps sin 9i)] (4 6) 

(Oy = \—9x cos 9 1 sm xpk 4- 9y cos 6j cos xpk + 9z sin 9i + 9i (Ps cos 6i + Pd sin 0/)] (4 7) 
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b)z = [9x sm 6j sm ipk - 6y sm 0/ cos tpk + 9z cos 9i + cos 9i + 9i{-Pd cos 9i + sm 9i)] 

(4 8) 

4.3 VELOCITY AT A POINT 

The absolute velocity vector F at a point ’p’ on the deformed beam can be written 
as 

F = Fff + rp + (w X f*p) (4 9) 

where Vh is the rigid body perturbational translation of the hub center Oh which 
IS given as 

Vh = ^1[Rx^xR + Ry^yR + Rz^zr] (4 10) 

Transforming all the unit vectors to 4K-system 

Vh = {VH)xex4k + {^H)yey4k + {yH)z&Zik (4 H) 

where 

{Vh)x= ^IRx {cos V’fc + sm ipk iPa cos 9 1 + Pd sin 9r ) } 

+ Q,lRy{cos tpkiPs cos 9j + Pd sm 9i) + sm '(j}k) 

+ QlRz{{Pp + PdCos9j - PsSin9i)} (4 12) 

{yH)y = 0.lRx{cosi)k{Ps <^os9i + Pd Sill 9 j — PpSin9i) — smi/)jfcCos^/} 

+ Q,lRy{cosipkCos9i + sixi‘ipkiPsCos9[ + PdSin9i - PpSm9i)} 

+ niRz{sin9i} (413) 
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(^''h)z — ^^Rxicosipki—PdCosOj + iSgSm 6 1 — /3pCos0/) + smtl^ksindj} 

+ ^lRx{cos'^k{—sm6i) + smtpk{—PcicosOi + / 3 s sin 0 / - ^ pcos ^/)} 

+ f ^ Z / 2 ^{ cos 0 /} (4 14 ) 

Tp = ex^M-adi{-psSm9i + pacosdr)+cosAsCOsAaUk 

+Vk sin As cos Aa — Wk sin A^ + cos Aj cos Aa cos(0g + (j)k) 

+u^(/i»fc cos As cos Aa sm{6a + ^k) - cos As cos Aa sin(0G + 4>k) 

—w'k cos As cos Aa<pk cos(0(j + ^k) — w'k cos As cos Aa Sin(0G! + (pk) 

+(f>k cos As sin Aa COS(0G + (^fc)} 

cos As COsAa COS(0G + (f>k) + w'k COS Aj COSAfl Sin(^G + 0fc) 
cos As cos Aa Sin(0G + ^fc) + v'k COS As cos AaPk cos{6g + (pk) 

+ sin As<Pk cos{6g + pk) - Pk cos As sin Aa sin(0G + Pk)}] 

+ey^i^lQ[—aOi sm6i + ei0/(— ^sSin0/ + P^cosOj — /3p cos 9 j) 

+620/ (—^S Sin 9 1 + cos 9i) - Uk sin As + Vk cos As 

sin As cos Aa COS(^G + Pk) + v'k sin As cos Aapk Sin(0G + Pk) 

—w'k sin As cos Aa Sin(0G + pk) - w'kPk sm As cos Aa COs(0G + Pk) 

-pk cos As sin(0G + Pk) + Pk sm As sm Aa cos(^g + Pk) } 

+C{-Wfc sm As cos Aa COS(0 g + Pk) + w'kPk sm As cos Aa Sm(^G + Pk) 

+v'k sm As cos Aa sin(0G + Pk) + v'k sm As cos AaPk cos{9g + 0k) 

-pk cos As sm(0G + Pk) - Pk sm As sm Aa sm(^G + Pk)}] 

+ez^JQ,[-a9f cos 0/ + ei0/(/3rf sm 9[ + Ps cos 0/ + Pp sm 9[) 

+ 620 / sm 9 1 + ;0s cos 9i) + cos A* sm Aa 
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sin As sm Aa + Wk cos Aq 

+77{i;^ sin Aa cos{da 4- 4>k) - v'^<pk sin(0G + <i>k) + sin Aa 
+ sin(0G + (/>*) + w'k4>k sin Aa cos{Bo + (j)k) + <Pk cos Aa cos{6g + 4>k)} 

+C{w',, sin Aa cos{6g + <j>k) - sin Aa<t)k sin(0G + (f)k) - v'^ sin Aa 
+ sin(0G + <f>k) - v'k sin Aa<j)k cos{6g + (j)k) - <t>k cos Aa sin(^G + (l>h)}] 

(4 15) 

and 

^X4k ^VAk ^^ik 

{uj xf) = oJx OJy Uz 

r* Ty Tz 

where rx,ry,rz and are the x,y and z components of fp and w re- 

spectively Substituting various quantities in Eq 4 9 from Eqs 4 2, 4 6-4 8 and 
4 11-4 16, the velocity at point p is obtained This can be written in symbolic form 
as 

V = ni[VxexAk + VyCy,, + VzSzJ (4 17) 

The detailed expression of the components of the velocity vector can be found 
in Fief [3] 
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Chapter 5 


EQUATIONS OF MOTION FOR 
ROTOR BLADE 


The coupled equations of motion of the rotor blade can be derived using Hamilton’s 
principle The mathematical form of Hamilton’s principle is stated as follows 

[^\6U - 5T - 5We)dt = 0 (5 1) 

Jh 

Where U is the strain energy, T is kinetic energy, Wg is the work done by the 
nonconservative external loads 

In this chapter, the expressions for the variation of kinetic energy and strain 
energy of the rotor blade are derived 


5.1 KINETIC ENERGY OF THE BLADE 


The total kinetic energy of the beam, T, is defined as 

T j j^pvyivoi 


(5 2) 
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Where V is the absolute velocity of an arbitrary point ’p’ on the blade cross 
section with respect to a reference frame The variation of kinetic energy, 5T can 
be written as 

= ' I dx (5 3) 


Substituting for the velocity V from Eq 4 17 and integrating (5T by parts with 
respect to time yields 




(ie). 


J p[Zu6u + Zy6v 4- Zjja6w + ZySv + Z'^Sw + Z^Sc^ildp dC dx (5 4) 


Where the terms Zv, Zyj,Zl, Z'^ and Z^ are the coefficients of du, 6v, 6w, 6v, 5w 
and 54> in the variation of kinetic energy expression 

Integration of the expression over cross-section yields 
6T = 52 f [Zu5u + Z^5v + Zy,6w + Zy5v + Z^^dw -I- Z^5^]dr} d( dx (5 5) 

I 

after eliminating the higher order terms, the expressions for Zu,Zy,Zju,Z'^,Z'^,Z<f, 
are obtained The detailed expressions of these terms can be found in Ref [3] 

The cross-section integrals are defined as 


m 

= IIa pdv d( 

mrjm 

= IIa PV dr] dC 

mCm 

= IIa pC dr] dC 

Imnn 

= / /a pC dvdC 


= /X4 PV^dpdC 

Irunc 

= IIa PV< dr] dC 


Where m is mass per unit length of the blade, mr)m and m(^m are first moments 
of cross-sectional mass per unit length, and IrUr,^ are the cross-sectional 

mass moments of inertia per unit length of the beam 
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5.2 TOTAL STRAIN ENERGY 


The formulation of strain energy obtained in this section essentially follows the 
procedure given in Refs [4] and [5] 


5.2.1 STRAIN ENERGY OF THE BEAM ELEMENT ( 
BLADE / FLEXBEAM / TORQUE TUBE ) 


The strain energy of a beam element is given by 



/ > 

T 

^ \ 

^xx 


^ XX 

* Txjj 

> < 

^XTf ^ 

7xC 

V > 


< 


dri dC dx 


(5 6) 


5. 2. 1.1 Explicit Strain-Displacement Relations 


The expression for non zero strain components written in terms of u,v,w and (j) are 
given in Refs [4] and [5] They are 

txx = Ux + \vx + 2^1 + + 04>i + 

+aTo{C‘ipr, - #c) “ [‘n<^os{6g + <!>)-( sm{dg + 

-[rj smiOg + 4>) + C cos{9g + (t>)]wxx + v{%r,,x “ + C{%c,x + '^ 0 %^) 

7x» 7 = 7i„ + - C(<^x + <^o) 

7xc = 7xc + ~ + <I>q) 

where 

4>q = {vxx COS do + vJxx Sin Og){~Vx sin Og + Wx cos da) 

The underlined term in Cxx represents the axial strain at the elastic axis These 
strain expressions can be simplified using the following assumptions 


• the transverse shear at the elastic axis is assumed to be zero 
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• The warping amplitude a is assumed to be equal to —(px 


The simplified strain components can be written as 
^xx =Ux + + \wl + 1 ( 77 ^ + - ^<l>xx - [roiC^T, - 

-[77 COs(0G +dk)-C Sm{9G + <I)k)]Vxx - [77 Sin(0G + ^k) + C COS(0g + <i)k)]Wxx 
'7x1} ~ (t/’jj d" C)*^i C*^0 

7xc = -(V'C + v)<px + V<f>Q 

(f>o = {vxx cos da + Wxx sm 9G)i-Vx sm 9g + Wx cos 9g) 

5. 1.1. 2 Stress-Strain Relations 

Assuming that the blade is made of isotropic material, the stress-strain relationship 
is given by the following equations 


f \ 

^xx 


0 

0 

i 


^XX 

^ ^XT) 1 

^ = 

0 G 0 

< 

7xij ^ 

% 

V./* 


I 

0 

0 


7x( 

V 


5.2.2 STRAIN ENERGY OF THE ELASTOMER 


The elastomer is placed between flexbeam and torque tube at a distance of Lgi from 
the root as shown m Fig 2 1 It is assumed to be placed between node on the 
flexbeam and node on the torque tube The elastomer is modelled as a combi- 
nation of three linear springs, one each in flap, lag and torsion modes respectively 

The strain energy of the elastomer Uei can be written as 


Uel = -KpiWi — WjY + — Vj)^ -b 

where Kp, Kl and are elastomer stiffnesses m flap, lag and torsion modes 
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5.2.3 STRAIN ENERGY OF THE CONTROL SPRING 


The control spring is located at a distance of Lc from the root with an offset Oc 
from the reference elastic axis It is attached at node on the torque tube The 
control spring is assumed to be linear having a stiffness Kc 

The deformation of the control spring due to deflection at the mode of the 
torque tube can be written as 


dc — + 0‘c^m 


The strain energy of the control spring Uc is given as 



(5 9 ) 


The total strain energy of the system can be written as U= C/^ + Uei + Uc 


5.3 TOTAL STRAIN ENERGY VARIATION 


5.3.1 VARIATION IN STRAIN ENERGY OF THE BEAM 
ELEMENT 

( BLADE / FLEXBEAM / TORQUE TUBE ) 


The variation in strain energy of a beam element is given by 


5Ub = Eol^ 

f \ 

T 

f N 

^ XX 


> < 

^XTJ ^ 


^7xC 


\ / 


dr] dx 


( 510 ) 
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5.3.2 VARIATION IN STRAIN ENERGY OF THE 
ELASTOMER 

The variation in strain energy of the elastomer SUei is given as 
6Uei = KF{'w^-WJ){6w^-6wJ)+KL{v^-VJ){5v^-5vJ)+K^{(j)^-(^J){6<f)^-5<|)J) (5 11) 

5.3.3 STRAIN ENERGY OF THE CONTROL SPRING 

The variation in strain energy of the control spring 6Uc is given as 

5Uc — Kc(Wjji + “h Q‘c^(pm) 


Total strain energy variation is (SU) = 6Ub + SUei + SUc 

The detail expressions of strains and their variations are available in Ref [3] 
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Chapter 6 


FORMULATION OF ELEMENT 
MATRICES ASSOCIATED 
WITH KINETIC AND STRAIN 
ENERGY VARIATION 

6.1 FINITE ELEMENT DISCRETIZATION 

The variational expressions associated with the kinetic and potential energy of the 
rotor blades are nonlinear The unknowns are the deformation functions u, v, w 
and <!> These are dependent on both space and time The spatial dependence is 
eliminated using Rayleigh-Ritz finite element formulation The blade is divided into 
subregions (finite elements) as shown in Fig 6 1, and the total dynamic potential is 
calculated for each subregion By applying Hamilton’s principle to each subregion, 
a discretized form of the equations of motion can be obtained In this development, 
each subregion is modelled by a straight beam type finite element These beam 
elements are located along the reference elastic axis (line of shear centres) of the 
blade 
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The discretized form of Hamilton’s principle is written as 



- 6T, - 5W^,)dt = 0 


(6 1 ) 


Where 


N= Total number of finite elements m the model 
5i7j= Variation of the strain energy of the element 
5Ti= Variation of the kinetic energy of the element 
SWei= Virtual work of the external loads on z*'* element 

The bearingless rotor is divided into three structural elements namely, (1) flexb- 
eam, (2) blade and (3) torque tube Each structural element is divided into finite 
number of beam elements Assume that the four unknown displacement functions 
of the beam finite element are expressed m the following form 


f N 

V 


’ {<l>vV 

0 

0 

0 


■ m ' 

w 

< 


0 


0 

0 


m 



0 

0 

{hV 

0 


{»} 

u 


0 

0 

0 

{<l>uV . 




(6 2 ) 


where are space dependent interpolation functions, and 

{V}, {W}, {$}, {U} are time dependent values of v,w, <j), u respectively 


{V} = 


✓ 


f N 

Wi 


f 
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<l>12 
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(t>2 
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\ ^ J 
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\ ^ J 



\ ✓ 


The nodal coordinates are shown m Fig 6 2 
The variation of the displacement function for the beam can be written as 


42 





6v 


ci 

O 

o 

o 

1 


■ {SV) ■ 

Sw 


o 

o 

o 


{W} 

S(f) 


0 0 0 


{«*) 

Su 

V J 


I 

J3 

o 

o 

o 


{«C/} 


In this development a cubic Hermite interpolation polynomial, {$c}) is used for 
the bending deflections (v,w) and a quadratic Lagrangian interpolation polynomial,{$g}, 
IS used for torsional rotation {(f>) and axial deflection {u) The mathematical expres- 
sions for these polynomials are given as 


1 1 1 1 


1 - 3^^ -b 2^^ 

ie{^-2e+e) 
se - 

hi-e+e) 



(6 4 ) 


{ } = { <^u } 


1 - 3^ + 2f 

4 ^ - 4^2 

K+a 



(6 5 ) 


where 

^ ~ ^k/^e 

Xk = spanwise (axial) coordinate of the beam element 
le = length of beam element 

For bending deformations, the nodal parameters are the displacements and slopes 
at both ends of the element Therefore, the resulting elements will have inter- 
element continuity for both displacement and slopes In addition, because of the 
cubic interpolation polynomial, bending strains vary linearly over the element length 
The quadratic Lagrangian interpolation polynomial is used for torsion(<^) and axial 
deformation(u) This polynomial has the capability of modelling a linear variation of 
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strains along the element length and therefore provides the same level of accuracy as 
the beam bending element The nodal parameters for these elements are chosen as 
the values of the displacements functions at the two end nodes and at the mid-point 
of the element 

The resulting beam element has 14 degree of freedom They are 4 m-plane (lag) 
bending degrees of freedom, 4 out-plane (flap) degrees of freedom, and 3 degrees of 
freedom each of torsion (<^), and axial deflection (u) The nodal degrees of freedom 
are shown in Fig 6 2 

6.2 ELEMENT MATRICES ASSOCIATED 
WITH KINETIC ENERGY VARIATION 

V 

The beam element matrices associated with the kinetic energy variation are derived 
by substituting the assumed expressions for the displacements functions (Eq 6 2) 
in the kinetic energy variation 6T (Eq 5 5) and carrying out the integration over 
the length of the beam element The resulting variation of the kinetic energy can 
be written in the form 


— — {59}^([M]i4x14{9} + [^^]iixu{Q} + [K‘^^]uxu{q} + {V’^}i4xi 


+[M^] 


14x3 


\ 

Rx 


r N 

Rx 


9^ 

* Ry 

> -f [M^]i4x3 * 

Ry 

► + [M®]i4x3 * 

Oy ^ 

Rz 

s ^ 


Rz 

K J 


9z 

V / 


-i-[M^]l4x3 S 


e. 

9,. 


9. 

\ j 


> + {V^^}l4xl + 


(6 6 ) 


Where {q} represents the vector of unknown nodal degrees of freedom 
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Detailed expressions for the various matrices defined m Eq 6 6 are given in Ref 


[3] 


6.3 ELEMENT MATRICES ASSOCIATED 
WITH STRAIN ENERGY VARIATION 


The elementary matrices associated with the strain energy can be obtained by taking 
variation of strain energy It can be written as 

W, = E„i={i?niAr®l{9} + {f®} («8) 

Where { q } represents the vector of unknown degree of freedom 
[K^\ IS the linear stiffness matrix {F^} is the nonlinear stiffness vector 
The detailed expressions for the [K^] and {F^} are found in Ref [3] 

6.4 ASSEMBLY OF ALL THE MATRICES 

After formulating all the elementary matrices for flexbeam, blade and torque tube 
the final assembly is performed The block diagram of global matrix is shown in 
Fig 6 3 

The contributions due to elastomer and control spring are added to the global 


matrix at the appropriate nodes 



(n-1) 1 e+ X k 










Al, A2, A3 and A4 together will form the junction element for flexbeam. 
Bl, B2, B3 and B4 together will form the junction element for torque tube. 
Cl, C2, C3 and C4 together will form the junction element for blade. 

Al, Bl and Cl are the independent dof. with junction, of flexbeam, blade 
and torque rube respectively. 

A2 and A3, B2 and B3 and C2 and C3 are the coupling terms of flexbam, 
blade and torque tube with junction dof. 

A4, B4 and C4 are the junction dof. of flexbem, blade and torque tube. 

Figure 6 3 Schematic diagram of global matnx 
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Chapter 7 


RESULTS AND DISCUSSION 


The first step in any aeroelastic response and stability analysis is the evaluation of 
natural frequencies and mode shapes of the rotor blade Using the inertial and struc- 
tural model developed m this study, a structural dynamic analysis was performed 
It may be noted that the inertial and structural operators given in equation (6 6) 
and (6 8) respectively, are non-linear Since the structural dynamic analysis requires 
only linear terms, all the non-hnear terms are neglected The corresponding linear 
equation for ith beam finite element can be written as 

(M1.{9}, + = 0 (71) 

Where [M]i represents the mass matrix of eth element, [K]^ represents the stiff- 
ness matrix of the ith element The detailed expressions for [M]i and [J^]i matrices 
can be found in Ref [3] The elementary matrices of flexbeam, blade and torque 
tube are assembled to get the global matrices as shown in Fig 6 3 

The contributions due to elastomer and control spring are added to the global 
matrix at the appropriate nodes Imposing the root boundary conditions, corre- 
sponding rows and columns are eliminated from the global matrices Then the 
resulting final global matrix equation can be written as 


[M]{q} + = 0 


(7 2) 
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Performing an eigen analysis, the natural frequencies of an undamped rotating 
blade in vacuum can be evaluated 


7.1 VALIDATION 

In order to validate the finite element blade model developed m this study, the 
results of the present analysis are compared (for certain specific cases) with those 
available m the literature The data shown m Table 7 1, correspond to a uniform, 
untwisted straight bearmgless rotor blade for which uncoupled natural frequencies 
have been evaluated m Ref [1] In this case the stiffness data for flexbeam and blade 
are same Torque tube is assumed to be rigid and massless The finite element code 
developed, m this study is suitably modified to get the natural frequencies for a 
bearmgless rotor blade with rigid and massless torque tube The terminology soft- 
in-plane blade configuration indicates that the first non-dimensional rotating lag 
frequency is less than 1 

Using the data given in Table 7 1, the coupled natural frequencies of the rotating 
bearmgless rotor blade are calculated The total number of rotor beam elements are 
20, of which 6 correspond to flexbeam and the remaining 14 correspond to blade 
The torque tube is attached to the flexbeam-blade junction at a distance(^i) of 0 Zl 
from the root The free end of the torque tube is at a distance(^ 2 ) of 0 from 
the root The elastomer is located at a distance(I/ei) of 0 1/ from the root In this 
case the control spring(lfc) has zero stiffness The results of the present analysis are 
shown in Tables 7 2 and 7 3 respectively for flap and lag modes 

The natural frequencies were evaluated for two different values of elastomer stiff- 
nesses It may be noted that even though the present formulation can treat coupled 
modes, for the assumed data the coupling between flap, lag and torsional modes is 
zero Prom Tables 7 2 and 7 3, it can be seen that the results of the present study 
are in excellent agreement with the results obtained in Ref [1] 
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7.2 effect OF MODELING OF TORQUE TUBE 

Influence of modeling the torque tube on the dynamic characteristics of the bear- 
ingless rotor blade has been analysed by treating the torque tube as either a rigid, 
massless beam (Model I) or as a flexible member having very high stiffness in all 
modes as compared to the stiffness of the flexbeam(Model II) 

The data pertaining to Model I is given m Table 7 1 Model II has the same 
data as Model I for flexbeam and blade But torque tube stiflPness values are as- 
sumed to be 500 times more than the properties of the blade in flap, lag and tor- 
sion In axial direction, torque tube stiffness is 50 times more than that of the 
blade Results are generated for two sets of different values of elastomer stiffness- 
es The non-dimensional stiffness values of the elastomer for the two cases are 
Kl = 30 1, Kp = 10 6 and Ki = 90 Z,Kp = 31 8 The natural frequencies of the 
rotor blade corresponding to the two torque tube models are shown in Table 7 4 

The results indicate that the natural frequencies are high for Model II as com- 
pared to Model I m flap and lag modes The influence of modelling the torque 
tube IS dominant in fundamental flap and lag modes The reason for this behaviour 
could be attributed to the strain energy stored m the torque tube for the Model II, 
whereas in Model I due to the assumption of rigid torque tube, there is no strain 
energy associated with it 

7.3 INFLUENCE OF CONTROL SYSTEM 
STIFFNESS 

For the study presented in Sec 7 1 and 7 2, the control system stiffness Kc has 
been set equal to zero To simulate a realistic bearmgless rotor blade, a control 
spring has been added to the rotor blade model A study has been performed 
to analyse the effect of control spring stiffness on the natural frequencies and the 
magnitude of coupling between flap-lag-torsion in various mode shapes The data 
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corresponding to this study are given in Table 7 5 The geometric parameters are 
6 = 0 3i, 6 = 0 01/, Lc = 0 02/ and = 0 05/ 

The influence of increasing control spring stiffness(i^c) on natural frequencies 
in first SIX modes are shown in Figs 7 1 to 7 6 Since control spring does not 
affect lag modes, natural frequencies m lag modes are not influenced by variation in 
control spring stiffness as can be seen in Figs 7 1 and 7 4 The other four natural 
frequencies exhibit an increase and reach asymptotic values with increase in Kc 

The effects of control spring on mode shapes are shown in Figs 7 7 to 7 12 
Lag mode shapes are not influenced by control spring as can be seen in Figs 7 7 
and 7 10 The effect of control spring is to highly couple the flap and torsional 
deformations, m flap modes as shown m Figs 7 8, 7 9 and 7 12 It is interesting to 
note that the participation of torsional motion substantially increases with increase 
m natural frequency m these modes 

The influence of control spring on torsion is shown m Fig 7 11 The effect 
of control spring seems to increase the torsional frequency from Ldn= 3 7095(with 
Kc = 0) to Un=- 4 4108(with Kc = 1000) with marginal influence on the mode 
shapes 


7.4 EFFECTS OF LOCATION OF ELASTOMER 
CONTROL SPRING AND FLEXBEAM- 
TORQUE TUBE JUNCTION 


The effect of varying the elastomer location, control spring location and the location 
of flexbeam-torque tube junction on the natural frequencies have been analysed for 
several combination of the geometric parameters The data used for this analysis is 
given m Table 7 5 the control spring stiffness is taken as Kc = 1000 The control 
spring location(Lc) is varied from 0 01/ to 0 03/ in steps of 0 01/ The elastomer 
location from the root(Le/) is varied from 0 05/ to 0 15/ in steps of 0 02/ The 
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junction, where flexbeam blade and torque tube meets, is located at a distance(^i) 
from the root and it is varied from 0 2ao 0 31 in steps of 0 02/ The offset of the 
control spring from the elastic axis is set at = 0 05/ The variation of natural 
frequencies are shown in 3D plots in Figs 7 13 to 7 18 

It IS evident from all these figures that the general behaviour of the variation 
of natural frequencies seems to be less influenced by the variation in location of 
control spring Lc The figures shows that all natural frequencies decrease as the 
elastomer moves away from the root This may be due to the increase m the length 
of the flexible element neai the root of the rotor blade The first two (1st lag, 
1st flap) natural frequencies exhibit a slightly increase and then decrease as the 
junction moves away from the root as shown m Figs 7 13 to 7 14 However the 
third and fourth natural frequencies(Figs 7 15 and 7 16) decrease as the location 
of the junction(^i) moves away from root On the otherhand, the sixth natural 
frequency (3rd flap) increases as the junction moves away from the root(Fig 7 18) 
The natural frequency in torsion increases as the junction moves away from the root 
For all these cases the numerical values are given in Appendix A 


7.5 EFFECTS OF TIP SWEEP AND ANHEDR- 
AL ANGLES WITH AND WITHOUT CON- 
TROL STIFFNESS 

The influence of tip sweep and anhedral angles on the natural frequencies and 
bendmg-torsion couplings of the rotor blade are analysed The bearingless rotor 
blade data is given in Table 7 5 The flexbeam, blade and torque tube are divided 
into 3, 7 and 4 elements elements respectively The swept tip portion of the blade 
is taken as one element having a length 0 11 The geometric parameters for all 
the cases are set at = 0 3/,^2 = 0 01/, Te/ = 0 05/ and I/c = 0 02/ The sweep 
and anhedral angles are varied from -30 to +30 in steps of 5 degree The coupled 
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natural fiequencies and mode shapes for two cases of control spring stiffness(/src == 0 
and Kc = 1000), have been evaluated In total there are 169+169 = 338 cases The 
variation of natural frequencies with sweep angles are shown in 3D plots from Fig 
7 19 to 7 24 

Figure 7 24 shows the variation of first lag frequency with sweep angles De- 
pending on the value of anhedral angle(Aa), the first natural frequency decreases 
monotomcally with sweep angle (for = + 30deg) or decreases and increases with 
increase in sweep angles (for Aq = -30deg) The first lag frequency varies from 
0 7533(for Ag = -30deg, Aa = 30deg) to 0 6396(for Ag = 30deg, Aq = 20deg) The 
control stiffness does not have any influence on lag frequencies A similar nature is 
also observed for second lag frequency, as shown in Fig 7 22 The second lag fre- 
quency varies from 3 3755(for Ag = -30deg, Aa = 30deg) to 3 1550(for Ag = 30deg, 
\a = 30deg) 

The variation of first flap frequency with sweep angle is shown in Fig 7 19 
Depending on the value of the anhedral angle the flap frequency increases monoton- 
ically(for Aa = +30deg) or decreases monotomcally (for Aa = -30deg) with increase 
in sweep angles The influence of sweep angle is more pronounced for large values of 
anhedral angles The first flap frequency varies from 1 1883 (for Ag = Aa = -30deg) 
to 1 11 54 (for Ag = -25deg, Ag = 30deg) without control spring When the control 
stiffness IS included the first flap frequency varies from 1 1942 (for A^ = Aa = 30deg) 
to 1 1201 (for Ag = 30deg, Aa = -25deg) A similar behavior can also seen in second 
and third flap frequencies as shown in Figs 7 21 and 7 24 

The torsional frequency seems to be less influenced by sweep and anhedral angles 
as shown m Fig 7 23 But the addition of control spring stiffness increases the 
torsional frequency 

It IS found from the results that anhedral angle alone introduces lag-torsion 
coupling and sweep angle alone introduces flap-torsion coupling 

The sweep angles introduce coupling between flap, lag and torsional deforma- 
tions in various modes To quantify the lag-torsion and flap-torsion couplings in 
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various modes, the value of torsional deformation at 0 751 ($o 75 ) is taken as refer- 
ence parameter, since this cross-section is treated as typical section for aerodynamic 
load calculations in a rotor blade The variation of $075 as a function of sweep 
angles is shown in Figs 7 25 to 7 30 

In first lag mode, for the case with Kc — 0, the lag-torsion coupling is significantly 
influenced by anhedral angle compared to sweep angle as shown m Fig 7 25a With 
the inclusion of control spring(/fc = 1000), both sweep and anhedral angle influences 
the lag-torsion coupling A similar observation can be seen in second lag mode shape 

It IS interesting to note that for the case with Kc = 0, the flap-torsion coupling 
in first flap is in the range -0 0760(for As = -30deg, Aq = Odeg) to 0 0760(for A* = 
30deg, Aa = Odeg) With the inclusion of control spring(A’c = 1000), the coupling 
range is -0 7967(for Ag = -30deg, Aa = -30deg) to -0 5859 (for A* = 30deg, Aa = 
-20deg) In case of second flap with Kc = 0, the flap-torsion coupling is m the range 
-0 0923(for A, = -30deg, Aa = -lOdeg) to 0 0923(for Ag = 30deg, A^ = lOdeg) With 
the inclusion of control spring(A’c = 1000), the coupling range is 2 8224(for As = 
-30deg, Aa = -30deg) to 3 0150(for A^ = 30deg, Aa = -15deg) In case of third 
flap with Kc = 0, the flap-torsion coupling is in the range -0 0832(for A* = 30deg, 
Aa = -5deg) to 0 0831 (for As = -30deg, Aa = 5deg) With the inclusion of control 
spring(i<’c = 1000), the coupling range is 20 3640(for A, = -30deg, Aa = -30deg) to 
26 9634(for A, = 30deg, Aa = -15deg) 

From the results it can be seen that the control spring introduces a large value 
of pitch-flap coupling as compared to the effect introduced by tip sweep angles the 
coupling parameter increases significantly with the increase in flap frequency 

The numerical values of frequencies and coupling parameters corresponding to 
Kc = 0 and Kc = 1000 are given in Appendix B 
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Coupled modes 


Ref [1], Uncoupled modes 


Spring 

stiffness 

Mode 

1 

FE Analysis 

No of elements = 20 

Power - Series 

Method 

Rayleigh- Ritz 

Method 


LAG 

1 

0 731 




LAG 

2 

4 453 



o 

11 

FLAP 

1 

1 133 

1 133 

1 134 

Kf = 10 6 

FLAP 

2 

3 432 

3 413 

3 416 


FLAP 

3 

8 072 

7 960 

8152 


TORSION 

1 

3 263 




AXIAL 

1 

6 939 




LAG 

1 

0 731 




LAG 

2 

4 453 

I 


o 

II 

FLAP 

1 

1 136 

1 135 

1 138 

Kf = 31 8 

FLAP 

2 

3 434 

3 415 

3 422 


FLAP 

3 

8 190 

8 075 

8 449 


TORSION 

1 

3 263 




AXIAL 

1 

6 939 




Cl = 0 3/,C2 = 0 ll,Lei = 01/ and Kc = 0 


Table 7 2 Comparison of natural frequencies in flap 
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Coupled modes 


Ref [1], Uncoupled modes 


■ II 


Spring 

stiffness 

Mode 

1 

FE Analysis 

No of elements = 20 

Power - Senes 

Method 

Rayleigh-Ritz 

Method 


LAG 

1 

0 780 

0 781 



LAG 

2 

4 459 

4 489 

4 493 

Kl = Z0l 

FLAP 

1 

1 125 



II 

o 

FLAP 

2 

3 427 




FLAP 

3 



7 715 






3 263 




AXIAL 

1 

6 939 




LAG 

1 

0 795 

0 796 

0 815 


LAG 

2 

4 460 

4 491 

4 499 

Kl = 90Z 

FLAP 

1 

1 125 



o 

II 

FLAP 

2 

3 427 




FLAP 

3 

7 715 




TORSION 

1 

3 263 




AXIAL 

1 

6 939 




^1 = 0 Zl, ^2 = 0 ll, Lei = 0 11 and Kc = 0 


Table 7 3 Comparison of natural frequencies in lag 
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Coupled inodes 




1 

MODEL I 

MODEL II 


■■ 

1 

iV= 10 

iV = 14 


LAG 

1 

0 780 

0 905 


LAG 

2 

4 459 

4 620 

= 30 1 

FLAP 

D 

1 133 

1224 

Kp = 10 6 

FLAP 

B 

3 432 

3 617 


FLAP 

B 

8 072 

8107 


TORSION 

1 

3 263 

3143 


AXIAL 

1 

6 939 

6 939 


LAG 

1 

0 795 

0 911 


LAG 

2 

4 460 

4 628 

Kl = 90 3 

FLAP 

1 

1 136 

1224 

Kf = 31 8 

FLAP 

2 

3 434 

3 634 


FLAP 

3 

8190 

8 268 


TORSION 

1 

3 263 

3143 


AXIAL 

1 

6 939 

6 939 


MODEL I Torque tube is assumed as rigid and massless 
MODEL II Torque tube is treated as a very stiff member 
Location of elastomer Lei = 0 li 

= 0 3Z,6 = 0 ll, Lei = 0 011 and Kc = 0 
Table 7 4 Effect of modelling of torque tube on blade natural frequencies 
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Soft-m-plane blade data 



Flexbeam 

Blade 

Torque tube 


N = 3 

iV = 7 

iV = 4 


0 0004 

0 0004 

0 0004 


00 

00 

00 


00 

00 

00 

m 

10 

09 

75 

Ps 


00 

00 

Pd 

00 

00 

ilDlHIl 

Pp 

00 



01 

00 

00 

00 

GJfmQ.n^ 

0 002 

0 003 

0 825 

EA/mOPP 

30 0 

40 0 

50 0 

e\ 

00 

0.0 

00 

62 

00 

00 

00 

a 

00 

00 

00 

EACq/EA 

0 00021036 

0 00021036 

0 00021036 


0 00206 

0 645 

0 505 

EInJmQPP 

0 0001 

0 006 

15 


Kp Elastomer stiffness m flap, 217 78 

Ki^ Elastomer stiffness in lag, 21 78 

Elastomer stiffness in torsion, 0 
Control spring stiffness, 1000 
^1=0 3^,6 = 0 U,Lc = 0 02i,ac = 0 051 and L^i = 0 051 


Table 7 5 Input data for soft-in-plane blade with flexible torque tube 
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Figure 7 2 Variation of second frequency ( 1st flap ) with control spring stiffness 
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Figure 7 4 Variation of fourth frequency ( 2nd lag ) with control spring stiffness 
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non dimensional blade deformation non dimensional blaae aetormai 


0 6656 


radial location 


Figure 7.7 (a) First natural frequency ( 1st lag ) mode shape with Kc=0 



COni = 0.6656 


radial location 

Figure 7 7 (b) First natural frequency ( 1st lag ) mode shape with Kc=1000 
Figure 7 7 Effect of control spring on mode shape of first natural frequency 





non dimensional blade deformation „ non dimensional blade del 


COmo ^ 1 1304 


radial location 

figure 7 8 (a) Second natural frequency ( 1st flap ) mode shape with Kc=0 



radial location 

Figure 7 8 (b) Second natural frequency ( 1st flap ) mode shape with Kc=1000 
Figure 7 8 Effect of control spring on mode shape of second natural frequency 




non 


*^‘^0 0.2 0.4 0.6 0.8 1 

radial location 

Figure 7 9 (b) Third natural frequency ( 2nd flap ) mode shape with Kc=1000 

Figure 7 9 Effect of control spring on mode shape of third natural frequency 
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non dimensional blade deformation non dimensional blade deformation 



Figure 7 10 (a) Fourth natural frequency ( 2th lag ) mode shape with Kc=0 



Figure 7 10 (b) Fourth natural frequency ( 2th lag ) mode shape with Kc=1000 
Figure 7.10 Effect of control spring on mode shape of fourth natural frequency 
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Figure 7 11 (a) Fifth natural frequency ( 1st torsion ) mode shape with Kc=0 



Figure 711 (b) Fifth natural frequency { 1st torsion ) mode shape with Kc-1000 
Figure 7 11 Effect of control spring on mode shape of fifth natural frequency 
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non dimensional blade deformation non dimensional blade deformation 



Figure 7 12 (a) Sixth natural frequency ( 3rd flap ) mode shape with Kc=0 
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25 
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0 


0.2 
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0)n6 ~ 5.9619 


0.8 


Figure 7 12 (b) Sixth natural frequency ( 3rd flap ) mode shape with Kc-1000 
Figure 7 12 Effect of control spring on mode shape of sixth natural frequency 
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FlgU! 


L 

e7l3 (a) 


01 



U 


el 

Control spring location, Lc=0 01 



Figure 7 13 (b) Control spring location, Lc=0 02 



Figure 713 (c) Control spring location, Lc=0 03 


Figure 7 13 Variation of first natural frequency ( 1st lag ) with location of elastomer 
and flexbeam-torque tube junction, for various control spring locations 



Figiue 7 14 (a) Control spring location, Lc=0 01 



Figure 7 14 (b) Control spring location, Lc=0 02 



Figure 7 14 (c). Control spring location, Lc-0 03 

Figure 7 14 Variation of second natural frequency ( 1st flap ) with location of 
elastomer and flexbeam-torque tube junction, for vanous control spring locations 
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FiguiP 7 15 (a) Control spring location, Lc=0 01 



Figure 7 15 (b) Control spring location, Lc=0 02 



Figure 715 (c) Control spnng location, Lc=0 03 

Figure 7 15- Variation of third natural frequency ( 2nd flap ) with location of elas- 
tomer and fiexbeam-torque tube junction, for various control spring locations 
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Figure 7 16 (a) Control spring location, Lc=0 01 



Figure 7 16 (b) Control spring location, Lc=0 02 



Figure 7 16 (c) Control spring location, Lc-0 03 

Figure 7 16 Variation of fourth natural frequency { 2nd lag ) with location of 
elastomer and fiexbeam-torque tube junction, for various control spring locations 
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Figure 7 17 (a) Control spring location, Lc=0 01 


4.5 

35 

Figuie 717 (b) Control spring location, Lc=0 02 



45 

4.0 

3.5 


Figure 7 17 (c) Control spring location, Lc—0 03 

Figure 7 17 Vanetion of fifth natural frequency ( 1st torsion ) with location of 
elastomer and frexbeam-torque tube junction, for various control spring locations 



'3 


Figure 7 18 (a) Control spring location, Lc=0 01 



Figure 7 18' (b) Control spring location, Lc=Q 02 



Figure 7 18 (c) Control spring location, Lc=0 03 

Figure 7 18 Variation of sixth natural frequency { 3rd flap ) with location of elas- 
tomer and flexbeam-torque tube junction, for various control spring locations 
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Figure 7 19 (a) for Kc=0 



Figure 7 19 (b) for Kc=1000 

Figure 7 19 Variation of first natural frequency! 1st lag )as a function of tip sweep 
angle with and without control spring stiffness 
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Figure 7 20 (a) for Kc=0 



Figure 7 20 (b) for Kc=1000 

Figure 7 20 Variation of second natural frequency ( 1st flap ) as a function of tip 
sweep angle with and without control spring stiffness 
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Figure 7 21 (a) for Kc=0 



Figure 7 21 (b) for Kc=1000 

Figure 7 21 Variation of third natural frequency ( 2nd flap ) as a function of tip 
sweep angle with and without control spring stiffness 
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Figure 7 22 (a) for Kc=0 



Figure 7 22 (b) for Kc=1000 

Figure 7 22 Variation of fourth natural frequency! 2nd lag ) as a function of tip 
sweep angle with and without control spring stiffness 
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i3 
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-30 


Figure 7 23 (a) for Kc=0 




Figure 7 23 (b) for Kc=1000 

Figure 7.23 Variation of fifth natural frequency{ 1st torsion ) as a function of tip 
sweep angle with and without control spring stiffness 
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Figure 7 24 (a) for Kc=0 



Figure 7 24 (b) for Kc=1000 

Figure 7 24 Variation of sixth natural frequency ( 3rd flap ) as a function of tip 
sweep angle with and without control spring stiffness 
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Figure 7.26 (b) for Kc=1000 


Figure 7 26 Effect of bendmg-torsion coupling as a function of tip sweep angle with 
and without control stiffness second frequency( 1st flap ) 








Figure 7.29 (a) for Kc=0 



Figure 7 29 (b) for Kc=1000 

Figure 7.29 Effect of bending-torsion coupling as a function of tip sweep angle with 
and without control stiffness fifth frequency( 1st torsion ) 
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Figure 7.30: (a) for Kc=0 



Figure 7 30- (b) for Kc=1000 

Figure 7 30. Effect of bending-torsion coupling as a function of tip sweep angle with 
and without control stiffness sixth frequency( 3rd flap ) 
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Chapter 8 


CONCLUDING REMARKS 


In this thesis, a linite element model has been developed for the dynamic analysis 
of a bearmgless lotor blade having a general configuration This model includes all 
geometric complelexities of the hub/blade and tip sweep 

The equations of motion have been derived using Hamilton’s principle The gen- 
eral equations aie then specialised for structural dynamic analysis of the rotor blade 
The formulation was validated by comparing the results of the present analysis for 
straight bearmgless rotoi blade, with those available m literature 

It is observed that the control spring stiffness introduces bendmg-torsion cou- 
pling The coupling parameter significantly increases with increase m flap natural 
frequency. Anhcdral angle introduces lag-torsion coupling and sweep angle intro- 
duces flap-torsion coupling. The results indicate that the control spring stiffness 
dominate the bending-torsion coupling as compared to the effect due to tip sweep 
angles 
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Appendix A 

Natural frequencies for variation 
of location of elastomer, control 
spring and flexbeam-blade-torque 
tube junction 



Table A 1 Natural frequencies for Lc — O.OlZ 


Lfii vanes fiom 0 05Z to 0 151 in steps of 0 021 
vanes from 0.21 to 0 in steps of 0 02Z 


^ei 

6 

1 lag 

1 flap 

2 flap 

2 lag 

torsion 

3 flap 

0.05 

0 20 

0.710 

1.145 

2.631 

3 617 

4.252 

5 448 

0 05 

0 22 

0.711 

1 147 

2 533 

3 483 

4.279 

5 522 

0 05 

0 24 

0 708 

1 147 

2.443 

3 383 

4 312 

5 608 

0 05 

0 26 

0 699 

1.145 

2 366 

3 307 

4 350 

5 703 

0 05 

0 28 

0 685 

1140 

2 299 

3.249 

4 379 

5 809 

0 05 

0 30 

0 666 

1.133 

2.247 

3 206 

4 401 

5.932 

0 07 

0 20 

0.712 

1141 

2 505 

3 096 

3 973 

5 444 

0.07 

0.22 

0.714 

1.142 

2.404 

2.986 

3.989 

5.512 

0.07 

0.24 

0.710 

1.141 

2.318 

2.905 

4.006 

5.590 

0 07 

0.26 

0.702 

1.138 

2.244 

2.846 

4 015 

5 678 

0 07 

0 28 

0 687 

1.132 

^2.184 

2.803 

4.015 

5.778 

0 07 

0 30 

0 668 

1.123 

2.136 

2.772 

3.999 

5.899 

0 09 

0 20 

0 712 

1.138 

2.402 

2 647 

3 824 

5 443 

0 09 

0 22 

0 713 

1138 

2 302 

2 549 

3 834 

; 5.507 

0 09 

0 24 

0 709 

1136 

2 216 

2 482 

3 838 

5 582 

0 09 

0 26 

0 699 

1131 

2146 

2 435 

3 833 

5 666 

0.09 

0 28 

0 683 

1124 

2 089 

2 402 

3 815 

5 765 

0 09 

0 30 

0 661 

1114 

2.046 

2 381 

3 785 

5 884 
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oil 0.20 0 711 1 135 2 323 2 291 3 740 5 444 

on 0 22 0.710 1.134 2 222 2 194 3 743 5 506 

on 0 24 0.702 1.131 2.138 2.134 3 740 5 578 

on 0 26 0 689 1 125 2 070 2 095 3 725 5 661 

on 0.28 0.670 1.117 2.017 2 071 3 698 5.758 

0.11 0 30 0.646 1 105 1975 2 057 3 658 5 877 


013 

0.20 

0 707 

1.133 

013 

0.22 

0.702 

1132 

0 13 

0 24 

0 691 

1 128 

0 13 

0 26 

0 673 

1.121 

0.13 

0.28 

0 649 

nil 

0 13 

0 30 

0 619 

1098 


0.15 

0.20 

0.703 

1.132 

0.15 

0 22 

0.693 

1.129 

0.15 

0.24 

0.674 

1124 

015 

0 26 

0 648 

1 116 

015 

0 28 

0 617 

1 105 

015 

0 30 

0.580 

1091 


2.260 

2 025 

3.687 

5 447 

2160 

1915 

3 687 

5 507 

2 077 

1855 

3 677 

5 577 

2 011 

1822 

3 657 

5 658 

1959 

1805 

3 623 

5 754 

1920 

1.798 

3 577 

5 873 


2 211 

1.851 

3 655 

5.450 

2.110 

1.708 

3 651 

5 509 

2.028 

1638 

3 637 

5 578 

1963 

1.606 

3 611 

5 658 

1.913 

1595 

3 573 

5 753 

1875 

1.594 

3 523 

5 872 
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Table A 2: Natural frequencies for Lc = 0.02i 


Lei vanes from 0 Qbl to 0 1 Si in steps of 0 02/ 
vanes from 0 21 to 0.3/ in steps of 0 02/ 


^el 

6 

1 lag 

1 flap 

0 05 

0.20 

0.710 

1.146 

0.05 

0.22 

0 711 

1 148 

0 05 

0 24 

0 708 

1.148 

0 05 

0 26 

0 699 

1.146 

0.05 

0.28 

0 685 

1 142 

0.05 

0.30 

0 666 

1 135 

0.07 

0 20 

0 712 

1 143 

0 07 

0 22 

0.714 

1.144 

0 07 

0 24 

0 710 

1143 

0 07 

0 26 

0.702 

1.140 

0 07 

0.28 

0.687 

1.134 

0 07 

0.30 

0.668 

1126 


0.09 

0.20 

0.712 

1.141 

0 09 

0 22 

0 713 

1 141 

0 09 

0 24 

0 709 

1 139 

0.09 

0 26 

0 699 

1.134 

0.09 

0 28 

0.683 

1127 

0 09 

0 30 

0.661 

1.117 


2 flap 

2 lag 

torsion 

3 flap 

2 639 

3 617 

4.255 

5 487 

2.540 

3.483 

4 283 

5 560 

2 449 

3 383 

4 317 

5 645 

2 371 

3 307 

4 357 

5 738 

2 304 

3 249 

4 387 

5 842 

2 250 

3 206 

4 411 

5 962 


2 515 

3 096 

3 972 
✓ 

5 484 

2 412 

2 986 

3 990 

5.551 

2 325 

2 905 

4.008 

5 627 

2 249 

2.846 

4 018 

5 713 

2.188 

2.803 

4.019 

5 810 

2.139 

2.772 

4 005 

5 927 


2.411 

2.647 

3.819 

5 483 

2 309 

2.549 

3 831 

; 5 546 

2 222 

2 482 

3 836 

5 619 

2.151 

2 435 

3.833 

5 701 

2.093 

2 402 

3 816 

5.795 

2 049 

2 381 

3 786 

5 911 
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Table A 3: Natural frequencies for Lc = 0 03/ 


Lei vanes from 0.05/ to 0 15/ in steps of 0 02/ 
vanes from 0 2/ to 0 3/ in steps of 0 02/ 


Ij^l 

6 

1 lag 

1 flap 

0 05 

0 20 

0.710 

1.148 

0 05 

0 22 

0 711 

1 150 

0 05 

0 24 

0 708 

1150 

0 05 

0 26 

0 699 

1.149 

0.05 

0.28 

0.685 

1.144 

0 05 

0 30 

0.666 

1 137 

0 07 

0 20 

0 712 

1.146 

0.07 

0.22 

0 714 

1.147 

0 07 

0 24 

0 710 

1.146 

0 07 

0 26 

0 702 

1143 

0 07 

0 28 

0.687 

1137 

0.07 

0 30 

0.668 

1.129 


0.09 

0.20 

0.712 

1.144 

0 09 

0 22 

0 713 

1.144 

0 09 

0 24 

0 709 

1.142 

0 09 

0 26 

0 699 

1138 

0 09 

0 28 

0 683 

1.130 

0.09 

0 30 

0 661 

1.121 


2 flap 

2 lag 

torsion 

3 flap 

2 649 

3 617 

4 258 

5 526 

2 548 

3 483 

4 287 

5 598 

2 456 

3 383 

4 323 

5 682 

2 376 

3 307 

4 364 

5 774 

2 308 

3 249 

4.395 

5 876 

2 254 

3 206 

4 420 

5 993 


2.525 

3 096 

3 971 

5 524 

2 421 

2 986 

3 989 

5 589 

2 331 

2 905 

4 009 

5 665 

2 254 

2.846 

4 021 

5 748 

2.193 

2 803 

4.023 

5.842 

2143 

2.772 

4 009 

5 956 


2.420 

2.647 

3.813 ' 

5.523 

2 316 

2.549 

3 826 

■ 5 585 

2 227 

2 482 

3 833 

5 656 

2155 

2 435 

3 831 

5 735 

2 096 

2 402 

3.816 

5.826 

2 051 

2 381 

3 786 

5.938 
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0 11 

0 20 

0 711 

1 142 

2 336 

2 291 

3 723 

0 11 

0 22 

0.710 

1 142 

2 232 

2 194 

3 730 

0 11 

0 24 

0.702 

1.139 

2.146 

2 134 

3 729 

0 11 

0 26 

0 689 

1 133 

2.076 

2 095 

3 717 

0 11 

0 28 

0 670 

1.125 

2 020 

2 071 

3 692 

0 11 

0 30 

0 646 

1.113 

1978 

2 057 

3 653 

013 

0 20 

0 707 

1.141 

2 270 

2 025 

3 667 

0.13 

0.22 

0.702 

1.140 

2166 

1.915 

3 670 

013 

0 24 

0 691 

1136 

2.081 

1855 

3 663 

013 

0 26 

0 673 

1129 

2 013 

1822 

3 645 

0 13 

0 28 

0 649 

1119 

1959 

1805 

3 613 

0 13 

0 30 

0 619 

1 107 

1919 

1798 

3 569 

015 

0 20 

0 703 

1 140 

2 217 

1 851 

3 631 

0.15 

0.22 

0.693 

1.138 

2.113 

1.708 

3 631 

0.15 

0 24 

0 674 

1.133 

2 029 

1638 

3 620 

015 

0 26 

0 648 

1 125 

1962 

1606 

3 597 

0.15 

0 28 

0 617 

1.115 

1911 

1.595 

3 560 

0.15 

0 30 

0 580 

1.101 

1.872 

1594 

3 511 


I 


5 524 
5 584 
5 651 
5 728 
5 817 
5 929 

5 526 
5 584 
5 650 
5 725 
5 813 
5 924 

5 529 
5 586 
5 650 
5 723 
5.810 
5 921 
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Appendix B 

Natural frequencies and coupling 
measures 
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lahl(‘ B 1 Natural frequencies with tip sweep angles variation for Kc = 0 



Aa deg 

1 lag 

1 flap 

2 flap 

2 lag 

torsion 

3 flap 

-30 0 

-30 0 

0 706 

1 188 

2 355 

3 272 

3 703 

4 791 

-25 0 

-30 0 

0 695 

1 186 

2 344 

3 255 

3 704 

4 768 

-20 0 

-30 0 

0 686 

1 183 

2 333 

3 241 

3 704 

4 744 

-15 0 

-30 0 

0 678 

1 180 

2 321 

3 229 

3 704 

4 719 

-10 0 

-30 0 

0 672 

1 175 

2 309 

3 219 

3 704 

4 695 

-05 0 

-30 0 

0 667 

1 170 

2 296 

3 211 

3 705 

4 671 

00 0 

-30 0 

0 664 

1 164 

2 283 

3 206 

3 705 

4 647 

05 0 

-30 0 

0 663 

1 158 

2 270 

3 203 

3 705 

4 623 

10 0 

-30 0 

0 663 

1 151 

2 256 

3 202 

3 704 

4 600 

15 0 

-30 0 

0 664 

1 143 

2 243 

3 202 

3 704 

4 578 

20 0 

-30 0 

0 667 

1 134 

2 230 

3 205 

3 704 

4 557 

25.0 

-30 0 

0 671 

1 125 

2 217 

3 210 

3 704 

4 536 

30 0 

-30 0 

0 676 

1 116 

2 205 

3 216 

3 703 

4 516 


-30 0 

-25 0 

0 709 

1 176 

-25 0 

-25 0 

0 698 

1 174 

-20 0 

-25 0 

0 689 

1 171 

-15 0 

-25 0 

0 681 

1 168 

-10 0 

-25 0 

0 674 

1 164 

-05 0 

-25 0 

0 669 

1 159 

00 0 

-25 0 

0 665 

1 155 

05 0 

-25 0 

0 662 

1 149 

10 0 

-25 0 

0 661 

1143 

15 0 

-25 0 

0 661 

1137 

20 0 

-25 0 

0 663 

1 130 

25 0 

-25 0 

0 666 

1123 

30 0 

-25 0 

0 669 

1 115 


2 319 

3 283 

3 705 

4 720 

2 311 

3 264 

3 705 

4 703 

2 302 

3 247 

3 706 

4 686 

2 293 

3 233 

3 706 

4 668 

2 284 

3 222 

3 706 

4 649 

2 274 

3 213 

3 706 

4 631 

2 263 

3 206 

3 706 

4 612 

2 253 

3 202 

3 706 

4 594 

2 242 

3199 

3 706 

4 576 

2 231 

3199 

3 706 

4 558 

2 221 

3 200 

3 706 

4 541 

2 210 

3 203 

3 705 

4 523 

2 200 

3 207 

3 705 

4 507 
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-30 0 

-20 0 

0 712 

1 164 

-25 0 

-20 0 

0 701 

1 162 

-20 0 

-20 0 

0 691 

1 160 

-15 0 

-20 0 

0 683 

1 157 

-10 0 

-20 0 

0 676 

1 154 

-05 0 

-20 0 

0 670 

1 150 

00 0 

-20 0 

0 665 

1 146 

05 0 

- 20.0 

0 662 

1 142 

10 0 

-20 0 

0 660 

1 137 

15 0 

-20 0 

0 659 

1 132 

20 0 

-20 0 

0 660 

1.127 

25 0 

-20 0 

0 661 

1 121 

30 0 

-20 0 

0 663 

1 116 

-30 0 

-15 0 

0 716 

1 153 

-25 0 

-15 0 

0 704 

1152 

-20 0 

-15 0 

0 694 

1 150 

-15 0 

-15 0 

0 685 

1148 

-10 0 

-15 0 

0 677 

1145 

-05 0 

-15 0 

0 671 

1 142 

00 0 

-15 0 

0 665 

1 139 

05 0 

-15 0 

0 661 

1 136 

10 0 

-15 0 

0 659 

1 133 

15 0 

-15 0 

0 657 

1 129 

20 0 

-15 0 

0 656 

1 125 

25 0 

-15 0 

0 657 

1 121 

30 0 

-15 0 

0 658 

1117 


2 289 

3 293 

3 707 

4 662 

2 283 

3 272 

3 707 

4 650 

2 277 

3 253 

3 707 

4 638 

2 270 

3 238 

3 707 

4 626 

2 263 

3 225 

3 707 

4 612 

2 255 

3 214 

3 707 

4 599 

2 247 

3 206 

3 707 

4 585 

2 239 

3 200 

3 707 

4 571 

• 

2 231 

3 197 

3 707 

4 557 

2 222 

3195 

3 707 

4 543 

2 214 

3195 

3 707 

4 529 

2 206 

3197 

3 707 

4 515 

2 197 

3 200 

3 707 

4 501 

2 264 

3 303 

3 708 

4 615 

2 260 

3 280 

3 708 

4 608 

2 256 

3 259 

3 708 

4 600 

2 251 

3 242 

3 708 

4 592 

2 246 

3 228 

3 708 

4 583 

2 241 

3 216 

3 708 

4 574 

2 235 

3 206 

3 708 

4 564 

2 229 

3 199 

3 708 

4 554 

2 223 

3 194 

3 708 

4 544 

2 217 

3191 

3 708 

4 533 

2 210 

3 190 

3 708 

4 522 

2 204 

3 191 

3 708 

4 511 

2 197 

3193 

3 708 

4 500 


98 



-30 0 

-10 0 

0 720 

1 144 

2 243 

3 313 

3 709 

4 578 

-25 0 

-10 0 

0 708 

1 143 

2 241 

3 288 

3 709 

4 575 

-20 0 

-10 0 

0 697 

1 141 

2 239 

3 265 

3 709 

4 571 

-15 0 

-10 0 

0 687 

1 140 

2 236 

3 246 

3 709 

4 567 

-10 0 

-10 0 

0 679 

1 138 

2 233 

3 230 

3 709 

4 562 

-05 0 

-10 0 

0 671 

1 136 

2 230 

3 217 

3 709 

4 556 

00 0 

-10 0 

0 665 

1 134 

2 226 

3 206 

3 709 

4 550 

05 0 

-10 0 

0 661 

1 132 

2 222 

3198 

3 709 

4 543 

10 0 

-10 0 

0 657 

1 130 

2 218 

3 192 

3 709 

4 536 

15 0 

-10 0 

0 655 

1 128 

2 214 

3188 

3 709 

4 529 

20 0 

-10 0 

0 654 

1 125 

2 209 

3186 

3 709 

4 521 

25 0 

-10 0 

0 653 

1 123 

2 205 

3185 

3 709 

4 513 

30 0 

-10 0 

0 654 

1 120 

2 200 

3186 

3 709 

4 504 

-30 0 

-05 0 

0 724 

1 136 

2 227 

3 323 

3 709 

4 549 

-25 0 

-05 0 

0 711 

1 135 

2 227 

3 295 

3 709 

4 549 

-20 0 

-05 0 

0 699 

1 135 

2 226 

3 271 

3 709 

4 549 

-15 0 

-05 0 

0 689 

1 134 

2 225 

3 250 

3 709 

4 548 

-10 0 

-05 0 

0 680 

1 133 

2 224 

3 233 

3 709 

4 546 

-05 0 

-05 0 

0 672 

1 132 

2 223 

3 218 

3 709 

4 544 


00 0 -05 0 0 666 1 131 2 221 3 206 3 709 4 541 

05 0 -05 0 0 660 1 130 2 219 3 197 3 709 4 538 

10 0 -05 0 0 656 1 129 2 217 3190 3 709 4 534 

15 0 -05 0 0 653 1 128 2 214 3.185 3 709 4 530 

20 0 -05 0 0 651 1 127 2 212 3182 3 709 4 525 

25 0 -05 0 0 650 1 125 2 209 3180 3 709 4 519 

30 0 -05 0 0 650 1 124 2 206 3 180 3 709 4 513 
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-30 0 

two 

0 728 

1 129 

-25 0 

00 0 

0 715 

1 130 

-20 0 

00 0 

0 702 

1 130 

-15 0 

00 0 

0 691 

1 130 

-10 0 

00 0 

0 681 

1 130 

-05 0 

00 0 

0 673 

1 130 

00 0 

00 0 

0 666 

1 130 

05 0 

00 0 

0 660 

1 130 

10 0 

00 0 

0 655 

1 130 

15 0 

00 0 

0 651 

1 130 

20 0 

00 0 

0 649 

1 130 

25 0 

00 0 

0 647 

1130 

30 0 

00 0 

0 646 

1 129 

- 30.0 

05 0 

0 732 

1 124 

-25 0 

05 0 

0 718 

1 125 

-20 0 

05 0 

0 705 

1 127 

-15 0 

05 0 

0 693 

1 128 

-10 0 

05 0 

0 682 

1129 

-05 0 

05 0 

0 673 

1130 

00 0 

05 0 

0 666 

1131 

05 0 

05 0 

0 659 

1 132 

10 0 

05 0 

0 654 

1 133 

15 0 

05 0 

0 650 

1 134 

20 0 

05 0 

0 647 

1135 

25 0 

05 0 

0 645 

1135 

30 0 

05 0 

0 644 

1 136 


2 215 

3 333 

3 709 

4 528 

2 216 

3 303 

3 709 

4 531 

2 217 

3 277 

3 709 

4 534 

2 218 

3 254 

3 709 

4 536 

2 219 

3 235 

3 709 

4 537 

2 219 

3 219 

3 709 

4 538 

2 219 

3 206 

3 709 

4 539 

2 219 

3 196 

3 709 

4 538 

2 219 

3 188 

3 709 

4 537 

2 218 

3182 

3 709 

4 536 

2 217 

3 178 

3 709 

4 534 

2 216 

3176 

3 709 

4 531 

2 215 

3175 

3 709 

4 528 


2 206 

3 342 

3 709 

4 513 

2 209 

3 310 

3 709 

4 519 

2 212 

3 282 

3 709 

4 525 

2 214 

3 258 

3 709 

4 530 

2 217 

3 237 

3 709 

4 534 

2 219 

3 220 

3 709 

4 538 

2 221 

3 206 

3 709 

4 541 

2 223 

3195 

3 709 

4 544 

2 224 

3 186 

3 709 

4 546 

2 225 

3180 

3 709 

4 548 

2 226 

3175 

3 709 

4 549 

2 227 

3172 

3 709 

4 549 

2 227 

3170 

3 709 

4 549 


100 



-30 0 

10 0 

0 737 

1 120 

2 200 

3 350 

3 709 

4 504 

-25 0 

10 0 

0 721 

1 123 

2 205 

3 316 

3 709 

4 513 

-20 0 

10.0 

0 707 

1 125 

2 209 

3 286 

3 709 

4 521 

-15 0 

10 0 

0 695 

1 128 

2 214 

3 261 

3 709 

4 529 

-10 0 

10 0 

0 683 

1 130 

2 218 

3 239 

3 709 

4 536 

-05 0 

10 0 

0 674 

1 132 

2 222 

3 221 

3 709 

4 543 

00 0 

10 0 

0 665 

1 134 

2 226 

3 206 

3 709 

4 550 

05 0 

10 0 

0 659 

1 136 

2 230 

3 194 

3 709 

4 556 

10 0 

10 0 

0 653 

1 138 

2 233 

3185 

3 709 

4 562 

15 0 

10 0 

0 649 

1 140 

2 236 

3 177 

3 709 

4 567 

20 0 

10 0 

0 645 

1 141 

2 239 

3 172 

3 709 

4 571 

25 0 

10 0 

0 643 

1 143 

2 241 

3 168 

3 709 

4 575 

30.0 

10 0 

0 642 

1 144 

2 243 

3166 

3 709 

4 578 

-30 0 

15 0 

0 741 

1 117 

2 197 

3 358 

3 708 

4 500 

-25 0 

15 0 

0 725 

1 121 

2 204 

3 322 

3 708 

4 511 

-20 0 

15 0 

0 710 

1 125 

2 210 

3 290 

3 708 

4 522 

-15 0 

15 0 

0 696 

1 129 

2 217 

3 264 

3 708 

4 533 

-10 0 

15 0 

0 684 

1 133 

2 223 

3 241 

3 708 

4 544 

-05 0 

15 0 

0 674 

1 136 

2 229 

3 222 

3 708 

4 554 


00 0 

15 0 

0 665 

1 139 

2 235 

3 206 

3 708 

4 564 

05 0 

15 0 

0 658 

1 142 

2 241 

3194 

3 708 

4 574 

10 0 

15 0 

0 652 

1145 

2 246 

3183 

3 708 

4 583 

15 0 

15 0 

0 647 

1148 

2 251 

3175 

3 708 

4 592 

20 0 

15 0 

0 644 

1 150 

2 256 

3169 

3.708 

4 600 

25 0 

15 0 

0 642 

1152 

2 260 

3.165 

3 708 

4 608 

30 0 

15 0 

0 640 

1 153 

2 264 

3162 

3 708 

4 615 
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-30 0 

20 0 

0 745 

1 116 

2197 

3 365 

3 707 

4 501 

-25 0 

20 0 

0 728 

1 122 

2 206 

3 327 

3 707 

4 515 

-20 0 

20 0 

0 712 

1 127 

2 214 

3 294 

3 707 

4 529 

-15 0 

20 0 

0 697 

1 132 

2 222 

3 266 

3 707 

4 543 

-100 

20 0 

0 685 

1 137 

2 231 

3 243 

3 707 

4 557 

-05 0 

20 0 

0 674 

1 142 

2 239 

3 223 

3 707 

4 571 

00 0 

20 0 

0 665 

1 146 

2 247 

3 206 

3 707 

4 585 

05 0 

20 0 

0 657 

1 150 

2 255 

3193 

3 707 

4 599 

10 0 

20 0 

0 651 

1 154 

2 263 

3 182 

3 707 

4 612 

15 0 

20 0 

0 647 

1 157 

2 270 

3173 

3 707 

4 626 

20 0 

20 0 

0 643 

1 160 

2 277 

3167 

3 707 

4 638 

25 0 

20 0 

0 641 

1 162 

2 283 

3162 

3 707 

4 650 

30 0 

20 0 

0 640 

1 164 

2 289 

3159 

3 707 

4 662 

-30 0 

25 0 

0 749 

1 116 

2 200 

3 371 

3 705 

4 507 

-25 0 

25 0 

0 731 

1 123 

2 210 

3 331 

3 705 

4 523 

-20 0 

25 0 

0 714 

1 130 

2 221 

3 297 

3 706 

4 541 

-15 0 

25 0 

0 699 

1 137 

2 231 

3 268 

3 706 

4 558 

-10 0 

25 0 

0 686 

1 143 

2 242 

3 244 

3 706 

4 576 

-05 0 

25 0 

0 674 

1 149 

2 253 

3 223 

3 706 

4 594 

00 0 

25 0 

0 665 

1 155 

2 263 

3 206 

3 706 

4 612 

05 0 

25 0 

0 657 

1 159 

2 274 

3 192 

3 706 

4 631 

10 0 

25 0 

0 651 

1 164 

2 284 

3181 

3 706 

4 649 

15 0 

25 0 

0 646 

1 168 

2 293 

3172 

3 706 

4 668 

20 0 

25 0 

0 643 

1 171 

2 302 

3 165 

3 706 

4 686 

25 0 

25 0 

0 641 

1 173 

2 311 

3.160 

3 705 

4 703 

30 0 

25 0 

0 640 

1 175 

2 319 

3157 

3 705 

4 720 
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-30 0 

30 0 

0 753 

1 117 

2 205 

3 376 

3 703 

4 516 

-25 0 

30 0 

0 734 

1 126 

2 217 

3 335 

3 704 

4 536 

-20 0 

30 0 

0 716 

1 134 

2 230 

3 300 

3 704 

4 557 

-15 0 

30 0 

0 700 

1 143 

2 243 

3 270 

3 704 

4 578 

-10 0 

30 0 

0 686 

1 151 

2 256 

3 245 

3 704 

4 600 

-05 0 

30 0 

0 674 

1 158 

2 270 

3 224 

3 705 

4 623 

00 0 

30 0 

0 664 

1 164 

2 283 

3 206 

3 705 

4 647 

05 0 

30 0 

0 656 

1 170 

2 296 

3192 

3 705 

4 671 

10 0 

30 0 

0 650 

1 175 

2 309 

3180 

3 704 

4 695 

15 0 

30 0 

0 646 

1 180 

2 321 

3171 

3 704 

4 719 

20 0 

30 0 

0 643 

1 183 

2 333 

3164 

3 704 

4 744 

25 0 

30 0 

0 641 

1 186 

2 344 

3159 

3 704 

4 768 

30 0 

30 0 

0 640 

1 188 

2 355 

3155 

3 703 

4 791 
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Table B 2 Natuial fiequencies with tip sweep angles variation for Kc = 1000 


<ieg 

cif'g 

1 lag 

1 flap 

2 flap 

2 lag 

torsion 

3 flap 

-30 0 

-30 0 

0 706 

1 194 

2 388 

3 272 

4 588 

6 213 

-25 0 

-30 0 

0 695 

1 192 

2 377 

3 255 

4 572 

6 186 

-20 0 

-30 0 

0 686 

1 189 

2 366 

3 241 

4 556 

6 159 

-15 0 

-30 0 

0 678 

1 185 

2 353 

3 229 

4 539 

6 133 

-10 0 

-30 0 

0 672 

1 181 

2 341 

3 219 

4 523 

6107 

-05 0 

-30 0 

0 667 

1 176 

2 328 

3 211 

4 506 

6 082 

00 0 

-30 0 

0 664 

1 170 

2 314 

3 206 

4 489 

6 058 

05 0 

-30 0 

0 663 

1 163 

2 301 

3 203 

4 472 

6 035 

10 0 

-30 0 

0 663 

1 156 

2 287 

3 202 

4 456 

6 013 

15 0 

-30 0 

0 664 

1 148 

2 274 

3 202 

4 440 

5 993 

20 0 

-30 0 

0 667 

1 139 

2 261 

3 205 

4 424 

5 973 

25 0 

-30 0 

0 671 

1 130 

2 248 

3 210 

4 409 

5 955 

30 0 

-30 0 

0 676 

1 121 

2 235 

3 216 

4 394 

5 938 

-30 0 

-25 0 

0 709 

1181 

2 352 

3 283 

4 539 

6136 

-25 0 

-25 0 

0 698 

1179 

2 344 

3 264 

4 528 

6 118 

-20 0 

-25 0 

0 689 

1 176 

2 335 

3 247 

4 516 

6100 

-15 0 

-25 0 

0 681 

1 173 

2 325 

3 233 

4 503 

6 081 

-10 0 

-25 0 

0 674 

1 169 

2 315 

3 222 

4 491 

6 063 

-05 0 

-25 0 

0 669 

1 165 

2 305 

3 213 

4 478 

6 044 

00 0 

-25 0 

0 665 

1 160 

2 295 

3 206 

4 465 

6 026 

05 0 

-25 0 

0 662 

1 154 

2 284 

3 202 

4 451 

6 009 

10 0 

-25 0 

0 661 

1 148 

2 273 

3199 

4 438 

5 992 

15 0 

-25 0 

0 661 

1 142 

2 262 

3199 

4 425 

5 976 

20 0 

-25 0 

0 663 

1135 

2 251 

3 200 

4 412 

5 960 

25 0 

-25 0 

0 666 

1127 

2 241 

3 203 

4 399 

5 945 

30 0 

-25 0 

0 669 

1 120 

2 230 

3 207 

4 387 

5 931 
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-30 0 

-20 0 

0 712 

1 169 

2 322 

3 293 

4 499 

6 077 

-25 0 

-20 0 

0 701 

1 167 

2 315 

3 272 

4 491 

6 066 

-20 0 

-20 0 

0 691 

1 165 

2 309 

3 253 

4 483 

6 053 

-15 0 

-20 0 

0 683 

1 162 

2 302 

3 238 

4 474 

6 041 

-10 0 

-20 0 

0 676 

1 159 

2 294 

3 225 

4 464 

6 028 

-05 0 

-20 0 

0 670 

1 155 

2 287 

3 214 

4 455 

6 015 

00 0 

-20 0 

0 665 

1 151 

2 279 

3 206 

4 445 

6 002 

05 0 

-20 0 

0 662 

1 147 

2 270 

3 200 

4 435 

5 989 

10 0 

-20 0 

0 660 

1 142 

2 262 

3197 

4 425 

5 976 

15 0 

-20 0 

0 659 

1 137 

2 253 

3195 

4 414 

5 963 

20 0 

-20 0 

0 660 

1 132 

2 245 

3195 

4 404 

5 951 

25 0 

-20 0 

0 661 

1 126 

2 236 

3197 

4 393 

5 939 

30 0 

-20 0 

0 663 

1 120 

2 228 

3 200 

4 383 

5 927 

-30 0 

-15 0 

0 716 

1 158 

2 296 

3 303 

4 465 

6 032 

-25 0 

-15 0 

0 704 

1157 

2 292 

3 280 

4 461 

6 025 

-20 0 

-15 0 

0 694 

1 155 

2 288 

3 259 

4 455 

6 018 

-15 0 

-15 0 

0 685 

1 152 

2 283 

3 242 

4 450 

6 010 

-10 0 

-15 0 

0 677 

1 150 

2 278 

3 228 

4 443 

6 002 

-05 0 

-15 0 

0 671 

1 147 

2 272 

3 216 

4 437 

5 993 

00 0 

-15 0 

0 665 

1 144 

2 266 

3 206 

4 430 

5 984 

05 0 

-15 0 

0 661 

1 141 

2 260 

3199 

4 423 

5 974 

10 0 

-15 0 

0 659 

1 137 

2 254 

3 194 

4 415 

5 965 

15 0 

-15 0 

0 657 

1 134 

2 247 

3 191 

4 407 

5 955 

20 0 

-15 0 

0 656 

1 130 

2 241 

3190 

4 399 

5 946 

25 0 

-15 0 

0 657 

1 126 

2 234 

3191 

4 391 

5 936 

30 0 

-15 0 

0 658 

1122 

2 228 

3193 

4 383 

5 927 
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-30 0 

-10 0 

0 720 

1 149 

-25 0 

-10 0 

0 708 

1 148 

-20 0 

-10 0 

0 697 

1 146 

-15 0 

-10 0 

0 687 

1 145 

-10 0 

-10 0 

0 679 

1 143 

-05 0 

-10 0 

0 671 

1 141 

00 0 

-10 0 

0 665 

1 139 

05 0 

-10 0 

0 661 

1 137 

10 0 

-10 0 

0 657 

1 135 

15 0 

-10 0 

0 655 

1 132 

20 0 

-10 0 

0 654 

1 130 

25 0 

-10 0 

0 653 

1 127 

30 0 

-10 0 

0 654 

1 125 

-30 0 

-05 0 

0 724 

1 141 

-25 0 

-05 0 

0 711 

1 140 

-20 0 

-05 0 

0 699 

1 139 

-15 0 

-05 0 

0 689 

1 139 

-10 0 

-05 0 

0 680 

1 138 

-05 0 

-05 0 

0 672 

1 137 

00 0 

-05 0 

0 666 

1 136 

05 0 

-05 0 

0 660 

1 135 

10 0 

-05 0 

0 656 

1 134 

15 0 

-05 0 

0 653 

1133 

20 0 

-05 0 

0 651 

1 131 

25 0 

-05 0 

0 650 

1 130 

30 0 

-05 0 

0 650 

1129 


2 275 

3 313 

4 439 

5 998 

2 273 

3 288 

4 437 

5 995 

2 271 

3 265 

4 434 

5 991 

2 268 

3 246 

4 431 

5 987 

2 265 

3 230 

4 427 

5 982 

2 261 

3 217 

4 424 

5 977 

2 257 

3 206 

4 419 

5 971 

2 253 

3198 

4 414 

5 965 

2 249 

3192 

4 409 

5 959 

2 245 

3188 

4 404 

5 952 

2 240 

3 186 

4 398 

5 945 

2 235 

3 185 

4 392 

5 938 

2 230 

3186 

4 386 

5 931 


2 259 

3 323 

4 418 

5 972 

2 258 

3 295 

4 418 

5 973 

2 258 

3 271 

4 418 

5 972 

2 257 

3 250 

4 417 

5 971 

2 255 

3 233 

4 416 

5 969 

2 254 

3 218 

4 415 

5 967 

2 252 

3 206 

4 413 

5 964 

2 250 

3 197 

4 411 

5 961 

2 248 

3190 

4 408 

5 957 

2 245 

3185 

4 405 

5 953 

2 242 

3182 

4 401 

5 949 

2 239 

3180 

4 397 

5 944 

2 236 

3180 

4 393 

5 939 
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-30 0 

00 0 

0 728 

1 134 

-25 0 

00 0 

0 715 

1 134 

-20 0 

00 0 

0 702 

1 134 

-15 0 

00 0 

0 691 

1 135 

-10 0 

00 0 

0 681 

1 135 

-05 0 

00 0 

0 673 

1 135 

00 0 

00 0 

0 666 

1 135 

05 0 

00 0 

0 660 

1 135 

10 0 

00 0 

0 655 

1 135 

15 0 

00 0 

0 651 

1 135 

20 0 

00 0 

0 649 

1 135 

25 0 

00 0 

0 647 

1 134 

30 0 

00 0 

0 646 

1 134 

-30 0 

05 0 

0 732 

1 128 

-25 0 

05 0 

0 718 

1 130 

-20 0 

05 0 

0 705 

1 131 

-15 0 

05 0 

0 693 

1 132 

-10 0 

05 0 

0 682 

1 134 

-05 0 

05 0 

0 673 

1 135 

00 0 

05 0 

0 666 

1 136 

05 0 

05 0 

0 659 

1 137 

10 0 

05 0 

0 654 

1 138 

15 0 

05 0 

0 650 

1 139 

20 0 

05 0 

0 647 

1 140 

25 0 

05 0 

0 645 

1 140 

30 0 

05 0 

0 644 

1 141 


2 246 

3 333 

4 402 

5 954 

2 248 

3 303 

4 404 

5 957 

2 249 

3 277 

4 406 

5 959 

2 249 

3 254 

4 408 

5 960 

2 250 

3 235 

4 409 

5 961 

2 250 

3 219 

4 410 

5 962 

2 250 

3 206 

4 411 

5 962 

2 250 

3196 

4 411 

5 961 

2 250 

3188 

4 410 

5 960 

2 249 

3 182 

4 409 

5 959 

2 248 

3178 

4 408 

5 957 

2 247 

3176 

4 406 

5 954 

2 245 

3175 

4 404 

5 951 


2 237 

3 342 

4 391 

5 941 

2 240 

3 310 

4 395 

5 946 

2 243 

3 282 

4 399 

5 951 

2 246 

3 258 

4 403 

5 955 

2 248 

3 237 

4 407 

5 958 

2 250 

3 220 

4 410 

5 962 

2 252 

3 206 

4 413 

5 964 

2 254 

3195 

4 415 

5 966 

2 255 

3186 

4 417 

5 968 

2 256 

3180 

4 419 

5 969 

2 257 

3175 

4 419 

5 970 

2 257 

3172 

4.420 

5 970 

2 258 

3170 

4 420 

5 970 
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-30 0 

10 0 

0 737 

1 124 

2 231 

3 350 

4 384 

5 933 

-25 0 

10 0 

0 721 

1 127 

2 236 

3 316 

4 390 

5 940 

-20 0 

10 0 

0 707 

1 130 

2 241 

3 286 

4 397 

5 947 

-15 0 

10 0 

0 695 

1 132 

2 245 

3 261 

4 403 

5 954 

-10 0 

10 0 

0 683 

1 135 

2 249 

3 239 

4 409 

5 960 

-05 0 

10 0 

0 674 

1 137 

2 253 

3 221 

4 414 

5 966 

00 0 

10 0 

0 665 

1 139 

2 257 

3 206 

4 419 

5 971 

05 0 

10 0 

0 659 

1 141 

2 261 

3194 

4 424 

5 977 

10 0 

10 0 

0 653 

1 143 

2 264 

3185 

4 428 

5 981 

15 0 

10 0 

0 649 

1 145 

2 267 

3 177 

4 432 

5 986 

20 0 

10 0 

0 645 

1 147 

2 270 

3172 

4 436 

5 989 

25 0 

10 0 

0.643 

1 148 

2 272 

3168 

4 439 

5 992 

30 0 

10 0 

0 642 

1 149 

2 274 

3166 

4 441 

5 995 

-30 0 

15 0 

0 741 

1 122 

2 229 

3 358 

4 381 

5 929 

-25 0 

15 0 

0 725 

1 126 

2 235 

3 322 

4 389 

5 939 

-20 0 

15 0 

0 710 

1 130 

2 241 

3 290 

4 398 

5 948 

-15 0 

15 0 

0 696 

1 134 

2 248 

3 264 

4 406 

5 957 

-10 0 

15 0 

0 684 

1 137 

2 254 

3 241 

4 414 

5 966 

-05 0 

15 0 

0 674 

1 141 

2 260 

3 222 

4 422 

5 975 

00 0 

15 0 

0 665 

1 144 

2 266 

3 206 

4 430 

5 984 

05 0 

15 0 

0 658 

1 147 

2 272 

3194 

4 437 

5 992 

10 0 

15 0 

0 652 

1150 

2 277 

3183 

4 444 

6 001 

15 0 

15 0 

0 647 

1 153 

2 282 

3175 

4 451 

6 008 

20 0 

15 0 

0 644 

1 155 

2 287 

3169 

4 457 

6 016 

25 0 

15 0 

0 642 

1 157 

2 291 

3165 

4 463 

6 023 

30 0 

15 0 

0 640 

1159 

2 295 

3162 

4 468 

6 029 
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-30 0 

20 0 

0 745 

1 120 

-25 0 

20 0 

0 728 

1 126 

-20 0 

20 0 

0 712 

1 131 

-15 0 

20 0 

0 697 

1 137 

-10 0 

20 0 

0 685 

1 142 

-05 0 

20 0 

0 674 

1 147 

00 0 

20 0 

0 665 

1 151 

05 0 

20 0 

0 657 

1 155 

10 0 

20 0 

0 651 

1 159 

15.0 

20 0 

0 647 

1 162 

20 0 

20 0 

0 643 

1 165 

25 0 

20 0 

0 641 

1 167 

30 0 

20 0 

0 640 

1 169 

-30 0 

25 0 

0 749 

1 120 

-25 0 

25 0 

0 731 

1 127 

-20 0 

25 0 

0 714 

1 135 

-15 0 

25 0 

0 699 

1 141 

-10 0 

25 0 

0 686 

1 148 

-05 0 

25 0 

0 674 

1 154 

00 0 

25 0 

0 665 

1 160 

05 0 

25 0 

0 657 

1165 

10 0 

25 0 

0 651 

1169 

15 0 

25 0 

0 646 

1 173 

20 0 

25 0 

0 643 

1 177 

25 0 

25 0 

0 641 

1 179 

30 0 

25 0 

0 640 

1 181 
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2 229 

3 365 

4 381 

5 930 

2 237 

3 327 

4 392 

5 941 

2 245 

3 294 

4 402 

5 953 

2 254 

3 266 

4 413 

5 965 

2 262 

3 243 

4 424 

5 977 

2 270 

3 223 

4 434 

5 989 

2 279 

3 206 

4 445 

6 002 

2 286 

3 193 

4 455 

6 014 

2 294 

3 182 

4 465 

6 027 

2 301 

3 173 

4 475 

6 039 

2 308 

3 167 

4 484 

6 051 

2 315 

3 162 

4 493 

6 063 

2 320 

3 159 

4 501 

6 074 


2 231 

3 371 

4 385 

5 933 

2 241 

3 331 

4 398 

5 947 

2 252 

3 297 

4 411 

5 962 

2 263 

3 268 

4 424 

5 977 

2 273 

3 244 

4 437 

5 993 

2 284 

3 223 

4 451 

6 009 

2 295 

3 206 

4 465 

6 026 

2 305 

3 192 

4 478 

6 044 

2 315 

3 181 

4 492 

6 061 

2 325 

3 172 

4 505 

6 079 

2 334 

3165 

4 518 

6 097 

2 343 

3 160 

4 530 

6115 

2 351 

3 157 

4 542 

6133 



-30 0 

30 0 

0 753 

1 121 

2 236 

3 376 

4 392 

5 941 

-25 0 

30 0 

0 734 

1 130 

2 249 

3 335 

4 407 

5 957 

-20 0 

30 0 

0 716 

1 139 

2 261 

3 300 

4 423 

5 975 

-15 0 

30 0 

0 700 

1 147 

2 275 

3 270 

4 439 

5 994 

-10 0 

30 0 

0 686 

1 156 

2 288 

3 245 

4 455 

6 015 

-05 0 

30 0 

0 674 

1 163 

2 301 

3 224 

4 472 

6 036 

00 0 

30 0 

0 664 

1 170 

2 314 

3 206 

4 489 

6 058 

05 0 

30 0 

0 656 

1 176 

2 328 

3192 

4 506 

6 082 

10 0 

30 0 

0 650 

1 181 

2 340 

3180 

4 524 

6106 

15 0 

30 0 

0 646 

1 186 

2 353 

3 171 

4 541 

6 131 

20 0 

30 0 

0 643 

1 189 

2 365 

3164 

4 558 

6157 

25 0 

30 0 

0 641 

1 192 

2 376 

3159 

4 574 

6183 

30 0 

30 0 

0 640 

1 194 

2 387 

3155 

4 590 

6 209 
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Tabic B 3 Bcnding-torsion coupling with tip sweep angles variation for Kc = Q 


A, deg Aq, deg 1 lag 1 flap 2 flap 2 lag torsion 3 flap 


-30 0 
-25 0 
-20 0 
-15 0 
-10 0 
-05 0 
00 0 
05 0 
10 0 
15 0 
20 0 
25 0 
30 0 

-30 0 
-25 0 
-20 0 
-15 0 
-10 0 
-05 0 
00 0 
05 0 
10 0 
15 0 
20 0 
25 0 
30 0 


-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 
-30 0 

-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 
-25 0 


0 059 
0 063 
0 066 
0 069 
0 071 
0 072 
0 072 
0 072 
0 071 
0 069 
0 066 
0 063 
0 059 

0 053 
0 056 
0 059 
0 061 
0 063 
0 064 
0 064 
0 064 
0 063 
0 061 
0 059 
0 056 
0 053 


-0 064 
-0 057 
-0 047 
-0 037 
-0 025 
-0 013 
0 000 
0 013 
0 025 
0 036 
0 047 
0 055 
0 062 

-0 067 
-0 060 
-0 050 
-0 039 
-0 027 
-0 013 
0 000 
0 013 
0 026 
0 038 
0 049 
0 059 
0 066 


-0 088 
-0 078 
-0 065 
-0 050 
-0 034 
-0 017 
0 000 
0 017 
0 032 
0 047 
0 060 
0 070 
0 079 

-0 090 
-0 079 
-0 066 
-0 051 
-0 035 
-0 017 
0 000 
0 017 
0 033 
0 049 
0 062 
0 073 
0 082 


0 210 
0 223 
0 233 
0 240 
0 245 
0 247 
0 247 
0 243 
0 237 
0 228 
0 216 
0 202 
0186 

0190 
0 200 
0 208 
0 214 
0 217 
0 219 
0 218 
0 214 
0 208 
0199 
0189 
0 176 
0 161 


0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 

0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 
0 947 


0 056 
0 051 
0 043 
0 035 
0 024 
0 013 
0 000 
-0 013 
-0 026 
-0 039 
-0 051 
-0 062 
-0 072 

0 063 
0 056 
0 048 
0 038 
0 026 
0 014 
0 000 
-0 014 
-0 028 
-0 042 
-0 054 
-0 066 
-0 076 
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-30 0 

-20 0 

0 044 

-0 070 

-0 091 

0 162 

0 948 

-25 0 

-20 0 

0 047 

-0 062 

-0 080 

0170 

0 948 

- 20.0 

-20 0 

0 049 

-0 052 

-0 067 

0 176 
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rahh’ B } Bf’n<iing-t(viHiou coupling with tip sweep angles variation for Kc = 1000 
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-10 0 

0 010 

-0 611 

3 007 

25 0 

-10 0 

0 008 

-0 600 

3 012 

30 0 

-10 0 

0 006 

-0 590 

3 015 

-30 0 

-05 0 

0 012 

-0 704 

2 894 

-25 0 

-05 0 

0 012 

-0 699 

2 909 

-20 0 

-05 0 

0 011 

-0 693 

2 923 

-15 0 

-05 0 

0 011 

-0 685 

2 938 

-10 0 

-05 0 

0 010 

-0 676 

2 952 

-05 0 

-05 0 

0 009 

-0 666 

2 965 

00 0 

-05 0 

0 009 

-0 656 

2 976 

05 0 

-05 0 

0 008 

-0 645 

2 987 

10 0 

-05 0 

0 007 

-0 634 

2 995 

15 0 

-05 0 

0 006 

-0 624 

3 002 

20 0 

-05 0 

0 005 

-0 614 

3 007 

25 0 

-05 0 

0 004 

-0 604 

3 009 

30 0 

-05 0 

0 003 

-0 596 

3 008 


0 046 
0 042 
0 039 
0 035 
0 031 
0 027 
0 023 
0 018 
0 013 
0 007 
0 002 
-0 004 

-0 on 


0 923 
0 923 
0 923 
0 924 
0 924 
0 924 
0 924 
0 925 
0 925 
0 925 
0 926 
0 926 
0 926 


24 560 
24 745 

24 935 

25 130 
25 328 
25 528 
25 728 

25 928 

26 126 
26 320 
26 511 
26 695 
26 873 


0 023 
0 021 
0 020 
0 018 
0 016 
0 014 
0 012 
0 009 
0 007 
0 004 
0 001 
-0 002 
-0 005 


0 924 
0 924 
0 924 
0 924 
0 924 
0 925 
0 925 
0 925 
0 925 
0 925 
0 926 
0 926 
0 926 


25 411 
25 509 
25 610 
25 713 
25 816 

25 919 

26 022 
26 123 
26 222 
26 319 
26 413 
26 503 
26 590 
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-30 0 

00 0 

-0 000 

-0 693 

-25 0 

00 0 

0 000 

-0 690 

-20 0 

00 0 

0 000 

-0 685 

-15 0 

00 0 

0 000 

-0 678 

-10 0 

00 0 

0 000 

-0 671 

-05 0 

00 0 

0 000 

-0 663 

00 0 

00 0 

0 000 

-0 654 

05 0 

00 0 

0 000 

-0 645 

10 0 

00 0 

0 000 

-0 636 

15 0 

00 0 

0 000 

-0 627 

20 0 

00 0 

0 000 

-0 619 

25 0 

00 0 

0 000 

-0 611 

30 0 

00 0 

0 000 

-0 605 

-30 0 

05 0 

-0 012 

-0 684 

-25 0 

05 0 

-0 012 

-0 682 

-20 0 

05 0 

-0 011 

-0 679 

-15 0 

05 0 

-0 011 

-0 674 

-10 0 

05 0 

-0 010 

-0 669 

-05 0 

05 0 

-0 009 

-0 662 

00 0 

05 0 

-0 009 

-0 656 

05 0 

05 0 

-0 008 

-0 648 

10 0 

05 0 

-0 007 

-0 641 

15 0 

05 0 

-0 006 

-0 634 

20 0 

05 0 

-0 005 

-0 627 

25 0 

05 0 

-0 004 

-0 621 

30 0 

05 0 

-0 003 

-0 616 


2 911 

-0 000 

0 924 

26 084 

2 924 

-0 000 

0 924 

26 092 

2 936 

-0 000 

0 925 

26 100 

2 948 

-0 000 

0 925 

26 107 

2 960 

-0 000 

0 925 

26 113 

2 970 

0 000 

0 925 

26117 

2 979 

0 000 

0 925 

26 120 

2 987 

0 000 

0 925 

26 122 

2 993 

-0 000 

0 925 

26 122 

2 998 

-0 000 

0 925 

26 121 

3 000 

-0 000 

0 925 

26 118 

2 999 

-0 000 

0 926 

26 114 

2 996 

-0 000 

0 926 

26 109 


2 924 

-0 024 

0 925 

26 565 

2 934 

-0 022 

0 925 

26 482 

2 944 

-0 020 

0 925 

26 395 

2 953 

-0 018 

0 925 

26 305 

2 962 

-0 016 

0 925 

26 213 

2 970 

-0 014 

0 925 

26 118 

2 976 

-0 012 

0 925 

26 022 

2 982 

-0 009 

0 925 

25 924 

2 985 

-0 007 

0 925 

25 825 

2 987 

-0 004 

0 925 

25 726 

2 987 

-0 001 

0 925 

25 628 

2 984 

0 002 - 

0 925 

25 531 

2 979 

0 005 

0 925 

25 435 
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-30 0 

10 0 

-0 025 

-0 677 

-25 0 

10 0 

-0 024 

-0 677 

-20 0 

10 0 

-0 022 

-0 676 

-15 0 

10 0 

-0 021 

-0 673 

-10 0 

10 0 

-0 020 

-0 670 

-05 0 

10 0 

-0 018 

-0 665 

00 0 

10 0 

-0 017 

-0 661 

05 0 

10 0 

-0 016 

-0 655 

10 0 

10 0 

-0 014 

-0 650 

15 0 

10 0 

-0 012 

-0 644 

20 0 

10 0 

-0 010 

-0 639 

25 0 

10 0 

-0 008 

-0 634 

30 0 

10 0 

-0 006 

-0 631 


-30 0 

15 0 

-0 037 

-0 671 

-25 0 

15 0 

-0 035 

-0 674 

-20 0 

15 0 

-0 033 

-0 675 

-15 0 

15 0 

-0 031 

-0 675 

-10 0 

15 0 

-0 029 

-0 674 

-05 0 

15 0 

-0 027 

-0 672 

00 0 

15 0 

-0 025 

-0 669 

05 0 

15 0 

-0 023 

-0 665 

10 0 

15 0 

-0 020 

-0 662 

15 0 

15 0 

-0 018 

-0 658 

20 0 

15 0 

-0 015 

-0 654 

25 0 

15 0 

-0 012 

-0 651 

30 0 

15 0 

-0 009 

-0 648 


2 931 

-0 048 

0 925 

26 850 

2 938 

-0 044 

0 925 

26 674 

2 945 

-0 040 

0 925 

26 493 

2 952 

-0 036 

0 925 

26 307 

2 958 

-0 032 

0 925 

26116 

2 963 

-0 027 

0 925 

25 923 

2 967 

-0 023 

0 924 

25 728 

2 970 

-0 018 

0 924 

25 533 

2 972 

-0 013 

0 924 

25 338 

2 972 

-0 007 

0 924 

25 144 

2 970 

-0 002 

0 924 

24 953 

2 966 

0 004 

0 924 

24 766 

2 959 

0 010 

0 924 

24 584 


2 933 

- 0.072 

0 925 

26 940 

2 937 

-0 065 

0 925 

26 674 

2 941 

-0 059 

0 925 

26 399 

2 945 

-0 053 

0 924 

26117 

2 948 

-0 047 

0 924 

25 831 

2 951 

-0 040 

0 924 

25 542 

2 953 

-0 033 

0 924 

25 251 

2 954 

-0 026 

0 924 

24 961 

2 954 

-0 018 

0 924 

24 674 

2 953 

-0 010 

0 924 

24 391 

2 950 

-0 001 

0 924 

24 114 

2 945 

0 008 

0 923 

23 845 

2 938 

0 017 

0 923 

23 584 
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-30 0 

20 0 

-0 048 

-0 667 

2 930 

-0 094 

0 925 

26 846 

-25 0 

20 0 

-0 045 

-0 672 

2 931 

-0 086 

0 925 

26 490 

-20 0 

20 0 

-0 042 

-0 676 

2 932 

-0 078 

0 924 

26 123 

-15 0 

20 0 

-0 040 

-0 679 

2 933 

-0 069 

0 924 

25 749 

-10 0 

20 0 

-0 037 

-0 680 

2 934 

-0 061 

0 924 

25 370 

-05 0 

20 0 

-0 035 

-0 681 

2 934 

-0 052 

0 924 

24 989 

00 0 

20 0 

-0 032 

-0 680 

2 935 

-0 043 

0 923 

24 608 

05 0 

20 0 

-0 029 

-0 679 

2 935 

-0 033 

0 923 

24 231 

10 0 

20 0 

-0 026 

-0 678 

2 934 

-0 023 

0 923 

23 860 

15 0 

20 0 

-0 023 

-0 676 

2 933 

-0 012 

0 923 

23 498 

20 0 

20 0 

-0 020 

-0 673 

2 930 

-0 000 

0 922 

23 147 

25 0 

20 0 

-0 016 

-0 671 

2 925 

0 012 

0 922 

22 810 

30 0 

20 0 

-0 012 

-0 669 

2 919 

0 025 

0 922 

22 487 

-30 0 

25 0 

-0 057 

-0 664 

2 922 

-0116 

0 925 

26 581 

-25 0 

25 0 

-0 054 

-0 672 

2 919 

-0 105 

0 924 

26 138 

-20 0 

25 0 

-0 050 

-0 679 

2 918 

-0 095 

0 924 

25 683 

-15 0 

25 0 

-0 047 

-0 685 

2 916 

-0 084 

0 924 

25 220 

-10 0 

25 0 

-0 044 

-0 689 

2 915 

-0 073 

0 923 

24 754 

-05 0 

25 0 

-0 041 

-0 693 

2 915 

-0 062 

0 923 

24 288 

00 0 

25 0 

-0 038 

-0 695 

2 915 

-0 051 

0 922 

23 827 

05 0 

25 0 

-0 035 

-0 697 

2 915 

-0 039 

0 922 

23 375 

10 0 

25 0 

-0 031 

-0 697 

2 914 

-0 026 

0 922 

22 935 


15 0 25 0 -0 027 -0 697 2 913 -0 012 0 921 22 513 

20 0 25 0 -0 023 -0 696 2 912 0 003 0 921 22 110 

25 0 25 0 -0 019 -0 695 2 909 0 018 0 921 21 730 

30 0 25 0 -0 014 -0 694 2 905 0 034 0 920 21 375 
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tills I 


-30 0 

30 0 

-0 065 

-0 662 

-25 0 

30 0 

-0 061 

-0 674 

-20 0 

30 0 

-0 057 

-0 684 

-15 0 

30 0 

-0 053 

-0 693 

-10 0 

30 0 

-0 050 

-0 701 

-05 0 

30 0 

-0 047 

-0 708 

00 0 

30 0 

-0 043 

-0 713 

05 0 

30 0 

-0 040 

-0 717 

10 0 

30 0 

-0 036 

-0 720 

15 0 

30 0 

-0 031 

-0 722 

20 0 

30 0 

-0 026 

-0 722 

25 0 

30 0 

-0 021 

-0 722 

30 0 

30 0 

-0 016 

-0 721 


2 909 

-0 135 

0 924 

26 165 

2 904 

-0 123 

0 924 

25 639 

2 900 

-0 no 

0 923 

25 101 

2 897 

-0 097 

0 923 

24 556 

2 895 

-0 084 

0 922 

24 011 

2 895 

-0 071 

0 922 

23 473 

2 895 

-0 057 

0 921 

22 947 

2 896 

-0 043 

0 921 

22 439 

2 897 

-0 027 

0 920 

21 957 

2 899 

-0 on 

0 920 

21 506 

2 900 

0 007 

0 919 

21 091 

2 901 

0 025 

0 919 

20 715 

2 901 

0 044 

0 918 

20 382 
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